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THE  NORMAL  MODES  OF  CAVITY  ANTENNAS 


ABSTRAC  T 


This  dissertation  is  primarily  an  attempt  at  arriving  at  an  under- 
standing of  the  resonance  properties  of  small  cavity  antennas  by  studying 
their  normal  modes.  These  modes  are  similar  to  the  free  oscillations 
of  other  oscillatory  systems.  Each  one  has  a characteristic  configuration 
and  a complex  exponential  time  dependence.  It  is  assumed  that  there  is 
an  infinite,  discrete  set  of  them.  There  is  one  fundamental  distinction 
between  them  and  others  more  commonly  investigated;  namely,  that  they 
occupy  an  infinite  volume  and  that  the  damping  is  a result  of  radiation. 

One  result  of  this  is  that  the  customary  statements  of  existence,  orthogo- 
nality, and  completeness  cannot  be  applied.  These  questions  are  not  con- 
sidered in  the  present  work.  The  concept  of  Q is  discussed  and  it  is  de- 
fined in  terms  of  the  complex  frequency. 

The  normal  modes  are  actually  transient  oscillations,  and  one  section 
is  employed  for  the  development  and  discussion  of  some  pertinent  time- 
dependent  electromagnetic  theory.  An  equivalence  theorem  is  first  de- 
veloped. This  states  that  the  field  in  a source-free  region  is  equivalent 

AF  18(600)85  486-7  iv 


to  that  produced  by  a double  layer  of  surface  currents,  electric  currents 
J = H x n,  and  magnetic  currents  K = n x E,  where  n is  the  outward 
normal.  This  much  of  the  theorem  is  identical  with  the  customary  forms. 
In  addition  to  the  boundary  values  it  is  necessary  to  include  initial  values 
if  the  boundary  values  are  not  known  for  all  past  times.  Thus,  if  initial 
values  of  the  fields  are  known  at  t = ta,  and  boundary  values  for  t £ ta, 
then  the  surface  currents  plus  volume  impulse  currents  which  “fire”  at 
t = ta  generate  an  equivalent  field.  These  volume  impulse  currents  are 
of  two  kinds,  electric  and  magnetic,  and  these  have  strengths  D and  B, 
respectively.  The  termination  properties  of  the  equivalent  currents  are 
also  discussed. 

The  behavior  of  waves  with  a complex  exponential  time  dependence 
is  the  second  topic  of  section  II.  The  envelope  of  such  a wave  is  an  ex- 
ponential curve  which  increases  in  the  direction  of  propagation.  In  a 
waveguide  the  phase  velocity  is  greater  than  the  envelope  velocity,  so 
that  individual  waves  advance  with  respect  to  the  envelope  and  hence 
the  waves  themselves  increase  as  they  propagate.  An  example  is  given 
to  show  that  this  does  not  lead  to  paradoxical  situations.  It  is  further 
shown  that  these  fields  in  a waveguide  have  qua  si- steady- state  character- 
istics with  respect  to  cutoff  phenomena. 

The  third  section  of  the  dissertation  contains  a theoretical  treat- 
ment of  the  normal  modes  of  radiating  cavities.  An  arbitrary  cavity 
opening  into  a half-space  is  considered,  and  by  a variational  method  a 
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stationary  formula  for  the  complex  frequencies  (eigenvalues)  is 
obtained. 

In  the  cavity  region  the  field  is  expanded  as  a sum  of  cavity  modes. 
This  expansion  is  discussed  in  detail  in  an  appendix,  because  most 
previous  discussions  of  this  are  incorrect  in  that  the  expansions  used 
are  not  complete . 

The  theory  greatly  simplifies  when  the  cavity  is  restricted  to  be 
a section  of  waveguide,  since  ordinary  waveguide  theory  can  then  be 
used.  Further  simplifications  are  possible  when  it  is  assumed  that  the 
aperture  field  is  just  thelowest  order  waveguide  mode,  and  that  the 
cavity  has  a high  Q.  Several  interesting  analogies  to  circuits  appear 
when  these  assumptions  are  made. 

A series  of  calculations  was  made  for  the  first  normal  mode  of  a 
square  waveguide  cavity.  The  cavity  has  a side  dimension  a,  depth  Z , 
and  is  filled  with  a medium  of  dielectric  constant  e r.  It  is  specified 

j>  a 

by  the  four  parameters  — , e r,  — , and  Q,  where  X.  is  the  free-space 

a a. 

wavelength  at  resonance.  Curves  of  — and  Q vs.  for  various  values 

of  e , are  shown.  It  is  shown  that  there  is  a minimum  value  of  Q as- 
r 

a . . ^ ... 

sociated  with  every  value  of  — and  that  specific  values  of  — and  e r are 

necessary  to  obtain  it.  The  minimum  Q was  computed  and  is  shown  in 

a graph  as  a function  of  The  accompanying  parameters  — and  e r 

are  also  shown.  It  turns  out  that  the  minimum  Q graph  is  a straight 

line  on  log-log  paper;  the  following  formula  closely  represents  it  for 
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y < 0.35: 


mm 


0.424(-J) 


The  effect  of  a lossy  dielectric  may  be  found  by  using  a complex 
dielectric  constant.  The  relation  of  the  dielectric  loss  factor  to  the  Q 
and  radiation  efficiency  is  discussed. 

The  last  section  of  this  work  deals  with  some  experiments  that 
were  made  to  check  the  theory.  A brief  review  of  the  methods  of  cor- 
relating measured  and  calculated  resonance  frequency  and  Q is  first 

made.  Then  the  results  of  measurements  on  three  square  waveguide 

l 

cavity  antennas  are  given.  The  antennas  had  — values  of  7.7,  2.0,  and 
0.4,  with  corresponding  e r values  of  3.80,  3.88,  and  8.41.  The  first 
six  modes  of  the  first  one  were  measured,  four  modes  of  the  second 
were  measured,  and  only  the  first  mode  of  the  last  one  was  measured. 
The  measured  resonance  frequencies  agreed  very  well  with  the  calcu- 
lated ones;  i.  e.,  within  a few  per  cent.  The  values  of  Q did  not  agree 
so  closely,  but  the  agreement  was  generally  within  30%,  which  is 
reasonable.  Various  factors  which  might  affect  these  results  are 
discussed. 
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THE  NORMAL  MODES  OF  CAVITY  ANTENNAS 


I.  INTRODUCTION 

In  recent  years  small  cavity  antennas  have  become  important 
for  a variety  of  purposes.  A typical  antenna  is  shown  in  Fig.  1.  It  con- 
sists of  a deformation  in  a conducting  screen,  together  with  some 

exciting  mechanism.  The  aperture 
• is  on  the  order  of  a quarter  of  a 

Q 

wavelength  in  its  greatest  dimension. 
These  antennas  always  have  marked 
resonance  properties,  which  are  un- 

( *p 

desirable,  since  an  antenna  is  more 
useful  when  it  can  be  operated  over 
a broad  band  of  frequencies.  • The 
resonance  properties  are  therefore 
of  considerable  interest,  and  it  is 
the  purpose  of  this  dissertation  to  discuss  them. 

The  antenna  characteristics  which  are  of  interest  are  the  resonance 
frequency  and  the  Q.  The  Q,  or  quality  factor,  is  a measure  of  the 
sharpness  of  resonance,  and  it  can  be  directly  connected  to  the  operating 
range  of  the  antenna  in  terms  of  the  standing  wave  ratio  on  the  input 
transmission  line.  (See  section  IV  below.) 


Fig.  1.  Cavity  antenna. 
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The  general  definition  of  Q is  as  follows: 


Q = 2ir 


Average  energy  stored 
Energy  dissipated  per  cycle,  at  resonance 


(1.1) 


If  the  system  has  more  than  one  degree  of  freedom  there  is  a Q defined 
for  each  mode,  and  this  holds  for  a system  with  an  infinite  number  of 
degree  of  freedom,  as,  for  example,  a closed  cavity.  However,  when  the 
concept  is  extended  to  antennas  a serious  difficulty  appears.  The  energy 
in  an  electromagnetic  field  is  defined  in  terms  of  volume  integrals,  and 
they  are  infinite  for  an  antenna  because  E^  decreases  as  l/r^,  but  the 
volume  increases  as  r^.  Hence  Q as  defined  by  eq.  (l.l)  is  infinite  for 
an  antenna. 

This  difficulty  is  in  principle  resolved  by  stating  that  by  stored 
energy  we  do  not  mean  the  total  field  energy,  but  only  that  part  of  it  in 
the  “local”  field,  as  opposed  to  the  radiation  field.  But  this  statement 
does  not  clear  up  the  trouble  — we  must  now  decide  what  we  mean  by 
local  field.  This  must  be  done  very  carefully;  since,  in  general,  super- 
position does  not  apply  to  energy,  and  it  cannot  be  arbitrarily  broken 
down  into  components.  In  the  case  of  a waveguide  there  are  well-de- 
fined orthogonal  modes,  and  superposition  may  be  applied  to  the  mode 
energies.  Schwinger1  has  used  this  property  to  show  that  in  the  neigh- 
borhood of  a discontinuity  in  a waveguide  the  total  field  energy  is  the 
sum  of  two  parts — that  in  the  propagating  mode,  and  that  stored  in  the 
higher  modes.  He  showed  that  the  latter  component  by  itself  has  all  the 
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properties  which  are  associated  with  stored  energy  of  circuits. 

2 

Chu  examined  the  antenna  problem  by  means  of  orthogonal 
' modes.  He  considered  an  omnidirectional  antenna  and  expanded  the 

field  outside  a sphere  surrounding  the  antenna  in  spherical  waves. 

Each  wave  could  be  replaced  by  an  equivalent  circuit  which  had  a well- 
defined  Q.  In  this  way  Chu  obtained  relations  between  the  antenna  size, 
Q,  and  gain.  He  assumed  very  idealized  conditions  between  the  antenna 
and  the  surrounding  sphere,  so  that  the  theory  gives  limits  of  perform- 
ance of  actual  antennas. 

3 

Counter  defined  stored  energy  as  the  difference  between  the  total 
energy  and  a “flow”  energy,  which  he  defined  as  the  energy  in  that 
plane  wave  which  has  the  same  Poynting  vector  as  the  field  under  con- 
r sideration.  This  gives  a basis  for  calculating  a finite  Q for  an  antenna 

\ * from  eq.  (1.1).  However,  it  can  be  shown  that  this  definition  does  not 

agree  with  Schwinger’s  in  the  waveguide  case. 

In  this  paper  a different  method  of  attack  is  adopted;  Q will  be 
calculated  by  finding  the  normal  modes  of  oscillation  of  the  antenna.  It 
is,  of  course,  well-known  that  the  transients  give  an  indication  of  the  Q. 
For  example,  if  a system  “rings”  for  a long  time  after  the  source  is 
removed,  then  it  has  a high  Q.  A more  precise  statement  of  this  is 
now  made  by  a brief  discussion  of  a simple  RLC  circuit. 

Consider  a series  RLC  circuit,  as  in  Fig.  2.  The  energy  definition 
(1.1)  leads  directly  to  the  value. 
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u)rlj 

a-^r- 


(1-2) 


where 


cor  = = resonance  frequency.  (1.3) 


one  normal  mode  of  oscillation.  In  this 


If  the  terminals  in  Fig.  2 are  connected  the  resulting  circuit  will  have 

R L C 

i Wv /W5'  |( 

mode  the  voltage  and  current  have  a time 
dependence  of  the  form  where  Pn 


Fig.  2.  Elementary  series  circuit. 


is  a complex  number: 
Pn  = iJn  + ^n* 


(1.4) 


The  real  and  imaginary  portions  of  the  complex  frequency  pn  can  be 
found  in  terms  of  the  circuit  constants;  one  way  to  do  this  is  as  follows. 
The  assumed  exponential  time  dependence- allows  us  to  keep  the  concept 
of  impedance;  the  transient  (normal -mode)  impedance  is  identical  to 
the  steady-state  impedance,  with  a;  replaced  by  pn.  Now  consider 
terminals  across  the  capacitor,  as  in  Fig.  3.  A consideration  of  the 


I 


* In  this  paper  the  symbol  ] - 4J- 1 is  used  exclusively. 
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voltages  and  currents  shows  that  the  impedance  looking  to  the  left  at 


the  terminals  must  be  the  negative  of  that  looking  to  the  right.  Thus 


-r£-7i  = - (R  + jPnL)- 
JPn  C 


(1.5) 


When  (1.4)  is  inserted  and  real  and  imaginary  parts  separated,  one 
obtains 


e .A 

n 2L 


(1.6) 


2 _ 2 , 2 

wn  " “r  n* 


(1.7) 


Now  define  the  normal  mode  Q: 


Q = -t- 

n 2# 


(1.8) 


It  is  evident  that  this  is  approximately  equal  to  the  steady- state  Q 
(eq.  (12)).  It  is  only  approximate  because  the  resonance  and  normal- 
mode frequencies  are  slightly  different,  but  the  difference  is  less  than 
1%  for  Q > 3.5. 

The  normal  mode  Q may  be  directly  connected  with  the  stored 

energy.  When  Q is  high  the  field  decays  very  slowly,  and  we  may 

-2if  _t 

reasonably  speak  of  an  average  stored  energy  that  dies  out  as  e . 


Call  this  average  energy  U.  Then 


U ~ e 


(1-9) 


!U  = _ 2if  n 

Bt 
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This  compares  with  eq.  (1.1)  for  the  steady- state  Q. 

We  shall  subsequently  assume  that  the  radiating  cavity  has  normal 
modes  which  have  a similarly  defined  Q.  The  subject  of  normal  modes 
in  open  regions  is  not  altogether  new;  in  particular,  oscillations  about 
a sphere  and  a spheroid  have  been  discussed  for  a long  time.  These 
bodies  were  examined  first  because  their  surfaces  coincided  with 
separable  coordinate  surfaces.  The  sphere  problem  is  treated  in  de- 

4. 

tail  by  Stratton,  who  finds  expressions  for  the  fields  as  well  as  the 
first  few  eigenvalues.  Schelkunoff5  gives  a summary  of  the  results 

* 

obtained  by  several  writers  on  the  prolate  spheroid.  This  problem  is 
of  considerable  interest  because  a wire  is  the  limiting  case  of  a pro- 
late spheroid. 

Oscillations  on  other  shapes  have  apparently  been  considered 
only  to  a very  minor  extent.  Schelkunoff  also  discusses  modes  on  a 
very  thin  wire  bent  to  an  arbitrary  curve,  and  makes  some  calculations 
for  a straight  wire.  Apart  from  this  discussion,  the  author  can  find  no 
reference  to  normal  modes  of  more  complex  antenna  structures. 

The  present  investigation  uses  a variational  method  to  obtain  the 
solution.  This  is  a powerful  technique  which  quickly  yields  a good 
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PI 

• approximation.  It  is  an  outgrowth  of  the  work  done  by  Schwinger  on 

discontinuities  in  waveguides.  It  consists  essentially  of  obtaining  an 

% 

integral  equation  for  an  aperture  field  (or  obstacle  current)  and  then 
properly  manipulating  the  equation  to  obtain  a stationary  formula  for 
some  parameter  of  interest.  The  techniques  have  been  applied  to  a 
wide  variety  of  field  theory  problems,  such  as  impedance  calculations, 
scattering  and  diffraction  problems,  and  propagation  problems.  Some 

e,  7,  a 

typical  references  to  this  work  are  listed  in  the  references.  The 

present  application  is  similar  in  some  ways  to  the  theory  used  to 

9 

compute  the  propagation  constant  for  a slotted  waveguide. 

t 

I 
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II.  ELECTROMAGNETIC  TRANSIENTS 


In  this  section  we  discuss  several  aspects  of  time -dependent 
electromagnetic  field  theory  which  are  interesting  and  useful  in  con- 
nection with  the  normal  mode  theory  in  the  next  section.  First  the 
time  dependent  equivalence  theorem  is  established.  Then  the  behavior 
of  fields  with  a time  dependence  e^*  e~ ^ ^ is  discussed. 

1.  EQUIVALENCE  THEOREM 

The  equivalence  theorem  is  a frequently  used  theorem  for  scalar 
and  vector  fields.  It  states  that  the  field  interior  to  some  surface  S is 
equivalent  to  that  produced  by  certain  source  distributions  on  S.  The 
theorem  is  intimately  connected  with  uniqueness  theorems  which  state 
that  functions  which  satisfy  appropriate  equations  are  uniquely  specified 
inside  S if  their  values  on  S are  given.  For  steady-state  electromag- 
netic fields  the  equivalent  sources  are  electric  and  magnetic  currents, 

10 

with  strengths  Hx  n and  n x E . The  uniqueness  theorem  states  that 
either  tangential  E or  H specifies  the  interior  field,  so  that  with  ap- 
propriate boundary  conditions  either  set  of  currents  can  be  made  suf- 
ficient by  itself.  Most  classical  statements  of  this  result  include  surface 
charges  as  well  as  currents.11  This  is  a redundant  formulation,  since 
the  surface  charges  and  currents  are  connected  by  the  continuity  equation. 

The  time  dependent  theorem  follows  from  the  steady-state  case 
by  the  use  of  Fourier  integrals.  However,  this  method  obscures  some 
points,  and  it  is  of  interest  to  examine  this  case  independently.  The 
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general  uniqueness  theorem  for  a time -dependent  electromagnetic 
field17  states  that  the  distribution  of  E and  H throughout  V at  the 
initial  time  ta  plus  either  tangential  E or  H on  the  bounding  surface 
for  t£  ta  specifies  the  field  in  V for  t £ ta.  An  equivalence  theorem 
will  now  be  set  up  in  which  the  field  is  produced  by  electric  and  (or) 
magnetic  currents  on  S plus  a volume  distribution  of  impulse  currents 
which  “fire”  at  t = ta.  The  standard  statements  of  this  theorem 
(Larmor-Tedone  formulas)13  use  surface  charges  as  well  as  surface 
currents,  and  they  go  back  in  time  to  t = -so.  There  is  no  provision 
for  initial  values.  The  method  of  proof  used  here  is  substantially  the 
same  as  that  used  by  Professor  V.  H.  Rumsey  to  prove  the  theorepn 
for  the  case  of  sinusoidal  time  dependence.* 

Let  there  be  given  a region  V bounded  by  S containing  a linear 
isotropic  medium  (Fig.  4).  The  medium  constants  are  continuous 

functions  of  position;  it  is- assumed  that  any  sharp  boundaries  are 

n replaced  by  thin  layers  where  the  con- 
stants vary  rapidly  but  continuously. 
Inside  V is  an  electromagnetic  field 
generated  by  sources  J.  The  fields  sat- 
isfy the  Maxwell  equations: 


* Classroom  presentation. 
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(2.1) 


V x e = -\i 


3 H 

3 t 


* 


V x H =e  JJi  +cr  E + J. 

— ' 3 t — - 


(2.2) 


For  a Green’s  function  we  use  the  field  due  to  an  impulse  current 
JQ.  This  current  is  an  impulse  in  both  space  and  time;  it  is  defined  as 

Jc  = q a (p  - Q)S(t  - t0).  (2.3) 


The  8 -function  is  the  unit  impulse  function,  defined  by 


/: 


b 

a 


f(t)  8 (t  - tQ)  dt  = 1 ' ° 


f(tn)  if  a < tQ  < b 
0 - • otherwise. 


(2.4) 


The  symbol  8 (P  - Q)  stands  for  the  triple  product, 


S (P  - Q)  = § (xp  - xQ)  S (yp  - yQ)  8 (zp  - zQ).  (2.5) 

The  impulse  current  and  the  fields  it  would  produce  are  discussed  in 
Appendix  I. 

Let  G0  (PtQ,  Qt,  q)  and  F0  (Pt0,  Qt,  q)  be  the  electric  and  mag- 
netic fields,  respectively,  at  point  P at  time  tQ  due  to  an  impulse  cur- 
rent JQ  at  point  Q at  time  t.  By  the  law  of  superposition,  the  electric 
field  of  an  arbitrary  current  distribution  J can  be  written  as14 

wU  J (Qt)  Go  (PtQ,  Qt,  1)  dvQ  dt,  (2.6) 

where  1 is  a unit  impulse  current  parallel  to  J,  and  J is  the  magnitude 
* nationalized  MKS  units  are  used  throughout  this  paper. 


AF  18(600)85  486-7 


10 


of  JT.  This  may  also  be  written  as 

E(pto)  = Jv  dvQ  J2  G0  (PtQ,  Qt,  J)  dt. 


(2.7) 


Because  of  the  finite  velocity  of  propagation,  the  upper  limit  on  the 
time  integral  can  be  reduced  from  + co  to  (t0  - (r/c)),  where  r is  the 
distance  from  P to  Q. 

The  fields  G0  and  FQ*  are  the  response  of  the  system  (i.  e.,  the 
region  V)  to  the  elementary  source  JQ.  They  satisfy  the  equations 


3 F» 

V x go  = - (J.  — =° 

B t„ 


(2.8) 


V x Fn  = e 


3G0 

Tt7 


+ a G0  + q 8 (P  - Q)  8 (t  - tQ). 


(2-9) 


Inside  V the  constants  m e & are  the  same  for  the  Green’s  function  as 
for  the  fields  (E.H);  outside  V they  are  as  yet  arbitrary. 

Now  form  the  volume  integral  of  ^ * (E  * F0  - G0  * H). 


/y  V • (E  x F0)  dv  = fy  (F0  • V x E - E • V x F0)  dvp 


+ q 8 (P  - Q)  8 (t 


” *o) 


dvi 


(2.10) 


/v  V • (G0X  H)  dv  = /v{  - H-  A - ao  • [<  E + j]  1 dvp. 

(2.11) 

* The  arguments  of  these  and  similar 'functions  are  sometimes  omittedfor  convenience.  -In  any  such  case 
it  wfll  be  clear  what  the  argument  is. 
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The  integrations  are  performed  over  points  P.  When  they  are  sub- 


tracted and  the  divergence  theorem  is  used  on  the  left-hand  side,  one 
obtains 

f f 3h 

Js  lo  - 9ox  H)*  nds  = JVM - F0-^=-)dv 

+ L e (G0i£.-E-=2)dv+  J Gq  . Jdv-q  . E (Q t)  S (t  - t0). 

* * 3 1 3 tQ  * * 

(2.12) 

| 

Now  integrate  over  time,  from  t = ta  to  t = co  . The  time  ta  is  the 
initial  time,  while  tQ  is  a variable  time.  In  the  final  result  (e.g., 
eq.  (2.30)),  tQ  is  the  observer’s  time. 

/’ 00  /*  r r 

ta  dt  4 ® * £o  - H)  • nds  -J  dt  - Folf) 


r®  r bf  3 go  /*°°  r 
+ Jta  dt  iv  6 (2o  3#  - Zft;  ) iv  + Jta  <1*  iv  5°  • 1 dv  - q • E(Qt0). 

(2.13) 

The  time  dependence  of  GQ  and  FQ  can  occur  only  as  (t  - tQ),  hence 

3 F 3 F„ 

_o  _ _o 
3 tQ  3 t 

and  (2.14) 

3 9o  _ 3 9o 

3t0  3 t 


Interchange  the  order  of  integration  in  the  first  two  volume  integrals 


and  use  (2.14). 

The  result  is 

r j00  f 

00 

-GO  . 

q • E(Qt0)  = 

dv  + Jve  - ' So. 

dv  + j 

dt  / G, 
+ Jm  — 

V lt=‘a  V 

*a 

xa  v 

- J ^ dt  (E  x FQ  - Gc  x H)  . n ds. 

(2.15) 
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At  the  upper  limit,  t = so  , F0  and  G0  are  zero,  since 
F0(Pt0,  Qt,  q)  = G0(Pt0,  Qt,  q)  = 0 for  tQ  < (t  + J).  (2.16) 


Hence 


q • E(Qt0)  = f - B(Pta)  • F0(Pt0i  Qta,  q)  dvp 
+ fy  D(Pta)  • G0(Pt0,  Qta,  q)  dvp 

/CO 

dt  f G0  (PtQP  Qt,  q)  • J(Pt)  dvp 
*a 

/CO 

dt  f E(Pt)  x F0(Pt0,  Qt,  q)  ■ n dsp 

ta  •'s 

00 

+ f dt  f Gq(P  tQ,  Qt,  q)  X H(Pt)  . n dsp. 
^t,  JS 


(2.17) 


To  interpret  these  integrals  we  need  the  reciprocity  theorem. 

The  mathematical  statement  can  be  quickly  obtained  from  eq.  (2.15), 

> 

if  we  let  the  time  interval  be  (-  03  , + 00  ) instead  of  (ta,  00  ),  and  if  we 
let  (E,  H)  be  the  fields  at  (Pt)  due  to  an  impulse  current  m at  (Mtj): 


E(Pt)  = Gq  (Pt,  Mtj,  m) 


(2.18) 


H(Pt)  = F0  (Pt,  Mtp  m). 


(2-19) 


Then,  we  obtain  from  (2.15), 


q • G0(QtQ)  Mtj,  m)  - fdtf  Gc(Pt  , Qt,  q)  • m S (P  - M)  S (t  - tj)  dvp 

* m " i r 


'-oo  » y 


f -n  Fp(Pt,  Mtj,  m)  ■ F0(Pt0,  Qt,  q)  dvp 

J v t = -°° 


(Equation  continued  on  next  page.) 
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00 


t = -00 


+ f e G0(Pt,  Mtj , m)  • G0(pt0,  Qt;q) 

JV. 

/CD 

dt  / fe>(Pt,  Mtj,  m)  x F0(Pt0,  Qt.  q) 
•'C  L 


dvt 


- G0(Pt0,  Qt,  q)  x F0(Pt,  Mtj,  m)J-  n dsp. 


(2.20) 


The  second  factor  in  each  of  the  two  volume  integrals  on  the  right-hand 
side  is  zero  at  the  upper  limit  by  eq.  (2.16),  and  because  the  first  factors 
have  the  source  and  observer  times  interchanged,  they  are  zero  at 
t = - oo , The  surface  integral  must  be  a constant  because  the  left-hand 
side  is  independent  of  S.  To  show  that  the  constant  is  zero,  let  S re- 
cede to  infinity. 

By  eq.  (2.16)  the  zero  intervals  of  GQ(Pt,  Mtj,  m)  and 
F0(Pt0,  Qt,  q)  overlap,  because  the  distances  from  P to  points  M 
and  Q are  infinite.  It  follows  that  the  surface  integral  is  identically 
zero  for  all  t*  The  final  result,  the  reciprocity  theorem,  is 


q • G0(Qt0,  Mtj,  m)  = m-  G0(MtQ,  Qtp  q).  (2.21) 


The  method  used  to  establish  this  theorem  is  an  extension  to  the  time- 
dependent  case  of  the  method  used  by  Schelkunoff. 

An  alternative  form  of  the  theorem  is  obtained  as  follows.  Let 
Gl(Pt0,  Qt,  q)  and  Fi(Pt0,  Qt,  q)  be  the  electric  and  magnetic  fields 
of  an  impulse  magnetic  current  of  strength  q.  These  fields  satisfy 

* This  result  is  a special  case  of  a more’ general  theorem:  ‘the  surface  integral  Ej  x ^ X g^)«  n ds 

is  zero  when  all  sources  of  both  fields  are  on  the  same  side  of  S. 
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(2.22) 


3 Fi 

VxGx  = -M  -=±-  q S (P  - Q)S  (t  - tQ) 

d t0 

Vx£i  =e  1^1  +cr  Gj.  (2.23) 

3 tQ 

Now  form  V • {G0(Pt,  Mtp  m)  x Fj  (P tQ,  Qt,  q) 

- Gj(Pt0,  Qt,  q)  x F0(Pt,  Mtj,  m)} 


and  integrate.  After  manipulations  identical  to  those  above,  one  obtains 


-q  - F0(Qt0,  Mtj,  m)  = m • Gj(Mt0,  Qtj,  q).  (2.24) 


The  integrals  on  the  right-hand  side  of  (2.17)  are  now  examined 
separately.  With  the  help  of  eq.  (2.24)  the  first  term  becomes 

f - B(Pta)  • F0(Pt0,  Qta,  q)  dvp  = q f G^Q^,  Pt&,  B)  dvp. 

(2.25) 

This  is  just  the  electric  field  at  (QtQ)  due  to  magnetic  current  im- 
pulses which  “fire”  at  t = ta;  they  are  distributed  according  to  the 
magnetic  induction  B which  existed  then. 

Use  eq.  (2.21)  on  the  second  term  of  (2.17)  to  obtain 

f D(Pta)  • G0(Pt0,  Q^,  q)  dvp  = q • f G0(Qt0,  Pta,  D)  dvp.  (2.26) 
JY  JY 

This  is  the  electric  field  at  (QtQ)  due  to  electric  current  impulses 
which  “fire”  at  t = ta;  they  are  distributed  according  to  the  electric 
displacement  D which  existed  then. 
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The  surface  integrals  are  similarly  transformed: 


- f dt  f E(Pt)  x F0(Pt0,  Qt,  q)  • n dsp  = 
ta  Js 

00 

q * f dt  f G1(QVPt,(nxE))dsp 

• a.  " r* 


(2.27) 


G0(Pt0,  Qt,  q)  x H(Pt)  • n dsp  = 


n))  dsp. 

(2.28) 


These  terms  represent  the  electric  field  at  Qt0  due  to  magnetic  cur- 
rents K = n x E and  electric  currents  J =.H  * n .distributed  on  S,  for 

t ^ ta. 

The  last  volume  integral  is 


f dt  f G0(Pt0,  Qt,  q)  • J dvp  = q • f dt  f G0(QtQ,  Pt,  J)  dvp. 

' + " *ir  » i »tr 

(2.29) 


V 


This  is  the  contribution  to  the  field  at  (Qt0)  from  actual  currents 
within  V.  The  first  terms,  eqs.  (2.25)  through  (2.28),  represent  con- 
tributions from  sources  outside  V and  from  sources  existing  at  t < ta. 

Notice  that  the  source  q multiplies  all  the  source  integrals,  as 
well  as  the  electric  field  E(QtQ),  in  (2.17).  Since  q is  arbitrary  we 
may  drop  it  and  write 
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E(Qto)=f  dt  J G0(Qt0,  Pt,  J)  dvp  + J G^Qto,  pta,  B)  dvp 


*a  V 


00 

+ / G0(Qt0,  Pta,  D)  dvp  + f dt  f G0(Qt0,  Pt,  (HX  n))  dS] 

» * x * c* 


ta  "S 


00 

+ f dt  J Gj(Qt01  Pt,  (nx  E))dsp. 


(2.30) 


ta  S 


The  equivalence  theorem  symbolized  in  eq.  (2.30)  may  be  stated 
as  follows.  If  a source-free  linear  isotropic  region  V bounded  by  S 
is  excited  by  external  sources,  then  the  internal  fields  are  reproduced 
identically  by  surface  electric  and  magnetic  currents  distributed  on  S 
according  to  H x n and  n x E,  respectively.  The  effect  of  pre-existing 
fields  may  be  accounted  for  by  volume  electric  and  magnetic  impulse 
currents  distributed  according  to  D and  B,  respectively,  which  occur 
at  the  initial  time. 

If  the  sources  J are  outside  V,  the  first  term  of  (2.30)  is  zero. 

If  J exists  only  inside  V,  then  the  first  three  terms  alone  must  give 
the  field  inside  V,  and  the  surface  currents  therefore  give  zero  to  in- 
terior points.  Since  the  inside  and  outside  may  be  interchanged,  we 
may  infer  that  the  surface  currents  give  zero  on  the  source  side  of  S. 
This  can  also  be  seen  from  the  fact  that  the  surface  currents  have 
precisely  the  values  needed  to  terminate  the  field.  (At  a surface  cur- 
rent layer,  tangential  E is  discontinuous  by  an  amount  IC  and  tangential 
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H,  by  an  amount  J;  but  K = n x E and  J = H x n,  so  these  terminate 
the  field.)* 

Another  interesting  condition  arises  when  there  are  no  sources 
after  t = ta;  there  is  only  a transient  field  existing,  which  must  get 
weaker  in  time.  We  may  remove  the  surface  S to  infinity.  In  such  a 
case  the  sources  on  S cannot  contribute  to  any  finite  point  within  a 
finite  time,  so  the  volume  impulse  currents  generate  the  complete 
field.  Consider  the  field  immediately  after  t = ta-  The  currents  have 
disappeared,  and  the  field  essentially  has  the  prescribed  initial  value. 
Since  the  impulse  currents  give  zero  for  t < ta,  we  have  a field  which 
at  all  points  of  space  jumps  from  zero  to  D,  B at  time  ta,  when  there 
are  impulse  electric  and  magnetic  currents  of  strength  D,  B.  Hence 
we  postulate  that  volume  impulse  electric  currents  of  strength  J pro- 
duce a discontinuity  (in  time)  in  D equal  to  J,  at  the  time  the  impulse 
currents  occur.  Similarly,  volume  impulse  magnetic  currents  of 
strength  K produce  a discontinuity  in  B equal  to  K.  These  postulates 
are  analogous  to  those  concerning  the  discontinuities  in  space  produced 
by  a surface  distribution  of  current. 

We  now  can  make  a more  general  statement  concerning  termina- 
tion of  the  field.  The  equivalent  surface  and  volume  currents  may  be 
regarded  as  surface  currents  on  the  hyper-surface  bounding  the  (xyzt) 

13 

* This  termination  property  maybe  used  to  set  up  the  equivalence  theorem.  See  Baker  and  Copson, 
also  Schelkunoff,  S.  A.,' Some  Equivalence  Theorems  <of  Electromagnetics  and  Their  Application  to 
Radiation  Problems,  B.S.T.J. , Vol.  XV,  January  1936.  p.  92. 
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volume  of  interest.  These  currents  generate  a unique  field  in  the  four- 
dimensional region,  and  terminate  it  at  the  surface.  We  may  make  this 
volume  finite  by  closing  it  with  the  hyper-plane  t = t^.  The  appropriate 
currents  to  put  on  this  plane  are  the  volume  impulse  currents 
J = - D,  = - B.  This  leaves  us  with  a set  of  surface  currents 
which  produce  a field  inside  some  volume  in  (xyzt)  space.  Outside  this 
hyper-volume  the  fields  are  zero. 

This  result  is  a generalization  of  the  three-dimensional  picture 

1 3 

obtained  by  Larmor  and  Tedone.  Undoubtedly  it  could  have  been  ob- 
tained more  compactly  and  more  elegantly  by  starting  with  the  four- 
dimensional formulation  of  Maxwell’s  equations. 

So  far,  the  external  boundary  conditions  on  the  equivalent  currents 
have  not  been  specified.  It  is  evident  that  they  are  arbitrary  because 
the  currents  give  zero  outside  S and  any  boundaries  may  be  assumed 
there  without  affecting  the  internal  fields.  If  we  choose  a short  circuit 
on  S (n  x G = 0 on  S)  for  the  boundary  condition,  then  the  surface 
electric  currents  will  be  “shorted  out’’  and  only  the  magnetic  currents 
will  be  left,  radiating  against  an  electric  conductor.  This  may  also  be 
seen  by  examining  eq.  (2.17).  In  the  last  term,  G0  x n will  be  zero  on 
S,  and  thus  there  is  no  contribution  from  the  electric  currents.  Simi- 
larly, open-circuit  boundary  conditions  (nx  IT  = 0 on  S)  leave  only 
electric  currents  radiating  against  a “magnetic  conductor.”  These 

ideas  suffice  to  prove  the  uniqueness  theorem  mentioned  above , 

( j 

I 

i 
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that  either  tangential  E or  H on  S plus  initial  values  of  E and  H 

* 

throughout  V uniquely  specify  the  field. 

As  an  example  of  the  use  of  equivalent  currents,  consider  a 
cavity  antenna  radiating  through  an  infinite  conducting  ground  plane, 

as  in  Fig.  5.  Suppose  that  the 
• field  varies  sinusoidally  with  time 

Q 

until  t = tj,  at  which  time  the  aper- 
ture is  suddenly  shorted  with  an 
electric  conductor.  What  will 
happen  to  the  fields  in  the  half- 
space ? 

To  answer  this  question,  use 
the  equivalence  theorem.  The 
volume  V is  the  half-space, 
bounded  by  the  plane  of  the  ground  plane.  It  is  assumed  that  the  initial 
time  ta  is  back  far  enough  so  that  the  contributions  from  the  volume 
impulse  currents  have  become  insignificant.  Hence  the  field  at  Q is 
equivalent  to  that  produced  by  electric  currents  J = H x n and  magnetic 
currents  K = n x E on  the  plane  z = 0.  The  unit  vector  n points  into 
the  cavity.  The  magnetic  currents  K are  zero  except  over  that  portion 
of  the  plane  where  the  aperture  is,  because  n x E = 0 on  the  ground 

* There  is  no  assurance  that  an  arbitrarily  specified  E or  H will  produce  a physically  possible  field. 

The  question  of  existence  is  a perplexing  one  and  becomes  important  when  the  boundaries  of  the 
region  include  sharp  edges. 


Fig.  5.  Cavity  antenna. 
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plane.  Thie  picture  is  greatly  simplified  when  a short  circuit  boundary 
condition  is*  adopted.  For  now,  as  shown  above,  we  are  left  with 
K = n x E?  fti  the  aperture,  radiating  against  an  infinite  conducting  plane 
screen.  But  now  ordinary  image  theory  can  be  used  and  it  follows  that, 
in  the  half-space,  the  fields  produced  by  K on  the  conducting  screen 
are  identical  to  the  fields  produced  by  2 K radiating  in  free  space. 

The  time  dependence  of  K is  of  the  form 


K ~ cos  cj  t U(t  j - t), 


(2.31) 


where  U(tj  - t)  is  the  unit  step  function, 


u(tl  - t)  = 


0 if  t > tj 


1 if  t < t. 


(2.32) 


To  obtain  the  field  at  Q we  can  find  the  field  due  to  an  infinitesimal 

source  at  P with  the  above  time  dependence,  and  integrate  over  the 

aperture.  The  fields  due  to  this  elementary  source  are  discussed  in 

Appendix  I.  It  is  shown  there  that  the  field  at  Q will  vary  sinusoidally 

with  time  until  t = ti  + — , at  which  time  it  will  abruptly  cease,  possibly 

c 

with  a spike,  depending  upon  the  phase  of  E when  the  shutter  is  closed. 
The  total  field  at  Q is  the  sum  of  contributions  from  all  points  in  the 
aperture.  Since  the  shutter  is  closed  instantaneously,  tj  is  the  same 


* Fbr  an  alternative  treatment  of  plane  boundaries,  leading  to  the  same  result,  see  Smythe,  W.  R., 
The  Double  Current  Sheet  in  Diffraction,  Phys.  Rev.,  72,  December,  1947.  p.  1066. 
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for  all  points.  Hence  in  a time  r = (Fig.  6)  the  field  at  Q will 


give  us  very  little  information  regarding  the  bandwidth  of  the  antenna. 

There  is  no  characteristic  exponential  time  decay  in  the  above 
problem  because  it  is  a very  special  case.  The  equivalent  currents 
were  on  a plane  and  they  were  all  stopped  instantaneously.  On  a more 
general  body  one  would  get  the  exponential  time  dependence.  For  ex- 
ample, if  the  terminals  at  the  base  of  the  feed  probe  were  suddenly 
shorted,  the  energy  in  the  cavity  would  be  expected  to  reverberate  and/ 
gradually  leak  out  the  aperture.  The  time  decay  in  this  case  would  be 
connected  with  the  bandwidth. 

2.  FIELDS  WITH  AN  EXPONENTIAL  TIME  DEPENDENCE 
An  oscillating  field  which  is  dying  out  in  time  may  have  a time 
dependence  of  the  form  e^**  53  Normal  modes  of  a lossy  network 


AF  18(600)85  486-7 


22 


and  normal  modes  in  a closed  cavity  with  a lossy  dielectric  have  this 


ii 


time  dependence.  In  the  next  chapter  it  is  assumed  that  the  normal 
modes  of  a radiating  cavity  also  have  this  time  dependence.  The  radi- 
ating cavity  differs  from  the ‘first  cases  because  there  is  propagation. 
We  now  examine  this  effect. 

The  one -dimensional  wave  equation  for  a field  component  is 


b 2e  _ J_  BZE  _ 
B z^  c^  B t^ 


(2.33) 


If  the  time  dependence  is  e^\  where  p = co  + , then  the  solution  to 

this  is  the  sum  of  two  waves,  traveling  in  the  positive  and  negative 
z-directions: 


E = A e 


j«(t-f)  -£(*-?) 


c + B e 


(2.34) 


Consider  the  first  term.  The  first  exponential  factor  represents  a 
wave  traveling  in  the  positive  z-direction.  The  second  factor  represents 
an  exponential  decay;  for  any  value  of  z,  E dies  out  as  e However, 

at  a given  value  of  t,  E builds  up  in  the  positive  z-direction,  according 
to  e L , The  field  therefore  appears  as  a sine  wave  which  builds  up 
exponentially  in  the  direction  of  propagation.  The  complete  wave,  in- 
cluding the  envelope,  moves  with  a phase  velocity  c.  The  second  term 
in  (2.34)  gives  a similar  wave  traveling  in  the  negative  z-direction. 

These  are  exactly  the  type  of  waves  that  would  be  obtained  on  an  infinite 
string  which  is  excited  by  a source  which  gets  weaker  in  time  exponentially. 
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These  fields  become  infinite  for  large  values  of  z,  and  difficulties 
with  respect  to  surface  integrals  at  infinity  might  be  expected.  It  is 
safest  to  say  that  the  transient  started  at  some  finite  time  in  the  past, 
and  because  of  the  finite  velocity  of  propagation,  it  must  be  zero  for 
z > Z.  The  maximum  value  of  the  transient  will  be  determined  by  the 
initial  conditions.  If  it  is  a spherical  wave  then  it  is  zero  for  r > R, 
and  the  maximum  value  will  decrease  as  l/R. 
a.  TRANSIENTS  IN  WAVEGUIDES 
The  fields  in  a lossless  waveguide  are  much  more  complicated 
than  the  simple  plane  wave  described  above.  A waveguide  is  a dis- 
persive region,  because  the  phase  velocity  depends  on  the  frequency. 
Hence,  one  would  expect  distortion  of,  say,  a transient  pulse,  as  it 
traveled  down  the  waveguide.  However,  normal  modes  exist  in  wave- 
guide regions.  For  example,  there  would  be  normal  modes  in  the 
region  shown  in  Fig.  7.  The  load  may  consist  of  a sheet  of  lossy 

Load 

Fig.  7;  Waveguide  region. 

material,  or  it  may  be  an  aperture  in  an  infinite  screen.  The  only  case 
for  -which  there  are  no  normal  modes  is  when  the  load  is  matched;  i.  e., 
produces  no  reflections. 

The  steady-state  solutions  for  waves  in  the  guide  are  formally 
correct  for  the  normal  modes,  but  the  frequency  & must  be  replaced 


Short 

Circuit 


Waveguide 
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by  the  complex  frequency  p,  and  some  new  interpretations  will  be 
necessary.  For  convenience  the  steady-state  waveguide  theory  is  now 
briefly  summarized.16  If  E(xyz) , and  H(xyz)  are  fields  in  the  wave- 
guide then  the  transverse  components  may  be  expressed  as  a sum  of 
waveguide  modes: 

X a y 2 — ^ 

E(xyz)  = 2 \ (z)  e (xy)  = 2 (a  e + a e V ) e (xy)  (2.35) 

— v — v — 

H(xyz)  = X lv  (z)  hy  (xy)  = X (£y  e V e V ) hy  (xy).  (2.36) 


The  vector  functions  ev  and  h^  have  x-  and  y-components  only;  they 
consist  of  two  separate  sets  apiece,  those  going  with  waveguide  modes 
which  have  no  Ez  (TE),  and  those  which  have  no  Hz  (TM).  The  mode 
functions  satisfy  the  two-dimensional  equation 


V2  ev  (xy)  + k 2 ev  (xy)  = 0 


(2.37) 


£ V X n = 0 


(2.38) 


on  the  waveguide  walls.  This  equation  with  boundary  condition  has 


solutions  only  for  values  of  kvl  which  are  discrete,  real,  and  positive. 
These  eigenvalues  are  the  cutoff  wave  numbers  for  the  various  modes. 
The  frequency  and  waveguide  propagation  constant  are  related  by 


2 _ ,2  . 2 
K V ^ y y • 


(2.39) 


, 2 2 

k = co  fj.  e 


(2.40) 
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The  propagation  constant  is  a complex  number: 


yv  = °-v  + j $v  > (2-4 1) 

where 

av  = attenuation  constant  for  they  ^ mode  (2.42) 

__  2 7T  j-i 

- — = phase  constant  for  the  v in  mode.  (2.43) 

The  magnetic  modes  are  similarly  defined,  with  an  appropriate  boundary 
condition. 

The  modes  are  orthogonal,  and  may  be  normalized  to  unity: 


/ * In  ds  = f hv  • h ds  = hvtJi  (2.44) 

JS  JS 

f ds  = °-  (2-45) 

J S 

The  integration  is  over  the  waveguide  cross  section.  The  following  re- 
lations may  be  shown. 

ev  (xy)  = hj,  (xy)x  z (2.46) 


+ 

+ 


J M 
J SJ  yfe 


TE 


= Z, 


= Z, 


TM 


= transverse  wave  impedance,  for  TE  modes 

_> 

= transverse  wave  impedance,  for  TM  modes 

(2.47) 

+ 

— . (2.48) 

by  by 
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If  the  load  is  at  z = 0,  the  load  impedance,  for  the  v ^ mode,  is  defined 

Vv  (0) 


as 


M°) 


The  functions  Vv  (z)  and  lv  (z)  obey  transmission  line 


equations  and  are  identical  to  voltage  and  current  on  a transmission 
line  of  characteristic  impedance  zv  and  propagation  constant  • 

Now,  if  E (xyz)  and  H (xyz)  are  normal -mode  fields,  we  use  the 
same  theory,  replacing  w by  the  complex  frequency,  p = + j £ . It 

is  of  interest  to  examine  the  propagation  constants  Jv  . From  eq.  (2.39), 


yv  = kv  - k 


V 


(2.49) 


Substitute  for  y and  k: 


- /3  2 + j2  a,/3  - k 2 - (oj  ^ ^ + j 2-p  co  ) /j.  e 


(2.50) 


where,  for  convenience,  the  subscripts  have  been  omitted.  Upon  separating 
real  and  imaginary  quantities,  one  obtains 

a /3  = (2.51) 


o,2-/32=k2-(&;  2 - £ 2 ) ^ e = (w2-w2+<f2)^e  , 

(2.52) 

where  co  Q = cutoff  frequency  of  v ^ mode.  This  last  substitution  for  k 
is  not  valid  unless  ^ and  e are  real.  Eqs.  (2.51)  and  (2.52)  may  be 
solved  for  cl  and  /3  ; the  results  are 


f e v/znr 

g 


(2.53) 


P = T 8» 


(2.54) 
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where 


co  r 2 1 2 


8=viil1"(^)  - Wl  V 1 ’ ( "C)  "W 


(2.55) 


and,  as  before, 


Q = 2^  * 


(2.56) 


The  parameter  g may  be* called  the  “normal-mode  factor”  for 

co 

the  waveguide.  Some  sample  curves  of  g vs.  are  shown  in 


Fig.  8.  If  g - 0,  then  a,  - 0 and 
Q.  - oo  . This  is  the  steady-state 
case,  and  g for  this  case  is  shown  on 
the  graph. 

.00 

When  (-£3=-)  < 1 the  normal  mode 
may  be  said  to  be  “above  cutoff”  (of 
the  v waveguide  mode)  and,  if  Q is 
high,  g essentially  coincides  with  the 
steady-state  curve.  The  approximation 


for  this  case  is 


Fig.  8.  Normal-mode  factor  vs.  frequency. 


' co  _ 2 ,a>c 

i - (— r)  • (-a r><  !• 


(2.57) 


When  (-jp)  > 1 the  normal  mode  is  “below  cutoff”,  and,  if  Q is  high, 
g is  approximately  given  by 


2Q 


b£>  > I- 


(2.58) 
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It  follows  then,  from  eqs.  (2.53)  and  (2.54),  that,  for  high  Q, 


a =8  ± 


if 

V'  - £?> 


/3  = t a)  jnJi  - Hr) 


CO  - ? 


,Cx) 


}>  (^r)<  1 


(2.59) 


± w ^Te  ^/(^-£)2  - l 


CJ  „ 7 

<-^r  - 1 


£^>i. 

• O)  7 


(2.60) 


When  (-^j-)  < 1 the  propagation  constant  fi  for  a normal -mode  wave  is 
essentially  the  same  as  that  for  a steady -state  wave  above  cutoff,  at  the 
frequency  " . Below  cutoff  the  attenuation  constant  a,  for  a normal  - 
mode  wave  is  essentially  the  same  as  that  for  a steady-state  wave  below 
cutoff,  at  the  frequency  &>  . In  addition  to  these  “quasi-steady-state” 
characteristics,  the  normal -mode  wave  also  has  attenuation  above 
cutoff,  and  propagation  below.  (Actually,  the  attenuation  is  negative, 
because  a and  fi  have  opposite  signs.  This  point  is  discussed  below.) 

These  effects  of  attenuation  and  propagation  are  relatively  minor. 
For  example,  the  wavelength  in  the  waveguide  is  given  by 


\ _ 2tt  _ 

g P " g ' 


(2.61) 


where  \Q  is  the  wavelength  of  a plane  wave  in  the  medium,  at  the 
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frequency  co  . Below  cutoff  g is  very  small,  and  X.g  is  very  large. 
Over  a limited  region  the  wave  nature  of  the  disturbance  would  not  be 
evident.  One  would  observe  a field  with  almost  uniform  phase,  just 
as  for  a steady-state  wave  which  is  below  cutoff. 

A wave  traveling  in  the  positive  z-direction  is  described  by 


E = ej("  t - P z)  t - a z 


(2.62) 


where  a is  a negative  number.  This  wave  increases  exponentially  in 
the  direction  of  propagation.  Upon  using  eqs.  (2.53)  and  (2.54),  eq.  (2.62) 


becomes 


E = e 


(2.63) 


where 


v,  =_ 


1 " JTTe  g 


(2.64) 


(2.65) 


The  phase  velocity  of  the  waves  is  v^,  and  we  may  call  V£  the  envelope 
velocity.  For  waveguide  modes  which  have  (-gj— ) < < 1 (far  above 
cutoff),  vj  and  v2  approach  the  velocity  of  light,  and  the  field  becomes 
a plane  wave,  as  described  by  eq.  (2.34).  Below  cutoff  g - 0 and  the 
phase  velocity  becomes  very  large,  whereas  the  envelope  velbcity 
becomes  small. 

In  general  the  envelope  velocity  V2  is  less  than  the  phase  velocity 
vj.  Hence,  if  any  one  particular  sinusoid  were  watched  it  would  appear 
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to  increase  as  it  moved  down  the  waveguide.  A similar  phenomenon 
can  be  observed  with  water  waves;  if  a group  of  waves  advances  into 
still  water,  the  individual  waves  will  be  seen  to  advance  through  the 
group. 

The  fact  that  these  waves  increase  in  the  direction  of  propagation 
may  appear  paradoxical,  since  it  is  opposite  to  the  steady-state  attenu- 
ation effect.  To  clear  this  up,  it  must  be  remembered  that  normal 
modes  are  characteristic  of  regions  where  there  are  reverberations, 
and  the  reflected  waves  will  take  care  of  any  apparently  paradoxical 
increases  of  field  strength.  A simple  example  will  illustrate  this.  Let 
the  region  be  a radiating  waveguide,  as  shown  in  Fig.  9.  Suppose  a 
waveguide  mode  which  is  below  cutoff  is  excited  at  the  aperture,  and 
the  circumstances  are  such  that  the  field  of  that  mode  is  essentially  in 
phase  along  the  length  of  the  guide.  The  wave  will  increase  toward  the 
rear  wall,  and  the  reflected  wave  will  in  turn  increase  toward  the 
aperture.  The  magnitude  of  the  transverse  electric  field  is  illustrated 
in  Fig.  10,  and  that  for  the  magnetic  field,  in  Fig.  11.  The  incident  and 
reflected  waves  are  denoted  by  arrows  and  are  shown  dashed.  The  sum 
of  the  two,  the  total  field,  is  the  solid  line.  The  transverse  electric 
field  must  be  zero  at  the  rear  wall  and  hence  the  reflected  electric 
wave  is  out  of  phase  with  the  incident  one.  The  reflected  magnetic 
wave  is  in  phase  with  the  incident  one. 


»* 


stow 
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The  total  field  in  the  guide  is  the  same  as  that  of  a steady-state 
waveguide  mode  which  is  cut  off  — it  decreases  away  from  the  aperture 
according  to  a hyperbolic  sine  or  cosine.  The  ratio  of  the  magnitudes 
of  E and  H is  controlled  by  the  aperture  impedance.  One  important 
distinction  must  be  noted  between  this  and  the  steady-state  case  — the 
aperture  fields  are  essentially  thos e reflected  from  the  rear  wall,  rather 
than  incident  fields,  and  this  introduces  a minus  sign  into  the  impedance 
expressions.  This  point  will  be  noted  again  in  the  next  section. 


j 
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III.  NORMAL  MODES  OF  A RADIATING  CAVITY 
1.  GENERAL  THEORY 

We  now  consider  electromagnetic  oscillations  in  the  open  region 

bounded  by  a cavity  and  an  infinite  con- 
ducting plane  screen,  as  in  Fig.  12.  The 
electric  field  satisfies  the  wave  equation 

VxVx  E + /ze  li-£  =0  (3.1) 

Bt2 


•Q 


and  the  boundary  condition 
n x e = 0 


(3.2) 


Fig.  12.  Source-free  radiating  cavity. 


on  the  ground  plane  and  cavity  wall.  In 
addition,  the  tangential  component  of  the 
electric  field  is  continuous  across  the 
aperture,  where  there  may  be  a dielectric  discontinuity.  We  also  as- 
sume that  the  original  source  of  the  field  was  inside  the  cavity,  so  that 
the  energy  is  propagated  outward. 

Now  assume  that  there  exists  a discrete  set  of  solutions  which 


j p_  t 

have  a time  dependence  of  the  form  e , where 


Pn  = £yn  + Kn- 

The  wave  equation  then  becomes 

VxV»E„-k|En=0.* 

jpht 

* From  here  on  the  time  dependence  e “ will  be  suppressed. 


(3.3) 
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where 


(3.5) 


ki  = e • 


The  possible  electric  fields  En(xyz)  are  eigenfunctions  and  the 
propagation  constants  kn  are  eigenvalues.  The  questions  of  existence, 
orthogonality,  and  completeness  of  these  modes  are  ignored  here.  The 
conventional  statements  are  hard  to  apply,  because  of  the  open  region 
and  the  fact  that  the  modes  increase  as  (e  ar)/r  as  r goes  to  infinity. 

For  purposes  of  computing  Qn  only  the  eigenvalues  are  necessary, 
and  a variational  method  will  be  employed  to  find  them.  In  the  preceding 
section  it  was  shown  how  the  tangential  electric  field  over  a closed  sur- 
face can  be  used  to  generate  the  fields  interior  to  the  surface.  This 

technique  is  used  to  set  up  the  variational  formulation  of  the  problem. 

1 2 

The  magnetic  field  in  the  the  two  regions,  Hn  in  the  half- space,  and  Hn 
in  the  cavity,  will  be  expressed  in  terms  of  z x En  in  the  aperture,  and 

f 1 f 2 

then  we  will  set  J En  * Hn  • z ds  = J En  x Hn  ■ z ds.  This  is  ob- 

viously  true  for  the  correct  field  En,  but  its  value  lies  in  the  fact  that 

values  of  kn  calculated  from  this  equation  are  stationary  with  respect 

to  first-order  variations  of  En  about  the  true  value.  (See  Appendix  III.) 

The  first  step  in  this  procedure  is  to  find  the  magnetic  field  in  the 

half-space.  In  section  II  it  was  shown  that  the  fields  in  the  half-space 

1 

are  the  same  as  those  produced  by  magnetic  currents  K = -2zx  E 

flowing  in  the  aperture,  and  radiating  into  free  space.  They  are  derivable 

ia 

from  an  electric  vector  potential  as  follows: 
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where 


E =V  x f 


1 e"jkr 

K -A ds. 

— 4irr 


(3.6) 


Then 


1 1 2 
H — |V  V • F + k F ] , 

~ jp  M “ 


(3.8) 


1 jpe  f 1 -jkr  , 1 _-jkr 

H =ii—  J K 5 V V • (K  ) ds.  (3.9) 

— 4i  A[ — r — r J 


The  subscript  n,  referring  to  the  nth  normal  mode,  has  been  omitted, 
and  it  is  understood  that  these  remarks  apply  to  any  of  the  modes.  In 
part  3 below  the  subscript  will  be  used  again. 

The  order  of  differentiation  and  integration  is  immaterial,  because 
the  integration  is  with  respect  to  points  P;  the  differentiation,  to  points 
Q.  We  may  now  write,  in  the  customary  matrix  notation: 

1 r 1 

H (Q,k)  = J T (Q,P,k)  K (P)  dsp,  (3.10) 


where 


K(P)  = z x E(P), 


(3.11) 


and  T is  a free-space  tensor  Green’s  function: 


1 -ipe  , VV  e ~jkr 

T (Q,P,k)  = (I  + — ) 


(3.12) 


* For  a detailed  d iscu'ssion  of  tensor  Green’s  functions  see  reference  3 It 
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I is  the  unit  matrix.  The  components  of  F have  simple  interpretations. 

1 

For  example,  T (Q,P)  is  (-2)  times  the  x-component  of  magnetic 
field  at  Q due  to  a y-oriented  unit  magnetic  current  element  at  P.  In 
common  with  other  Green’s  functions  it  has  certain  symmetry  properties; 
i.  e., 

1 1 

r (Q.p)  = r (p,q).  (3.13) 


The  tilde  indicates  transpose.  This  follows  immediately  from  the 
reciprocity  theorem,  eq.  (2.21),  which  can  be  rewritten,  for  this  appli- 
cation,* as 

qrx  (Q,M)»  =m  F j (M,Q)q.  (3.14) 


If  the  transpose  of  the  right-hand  side  is  taken,  (3.13)  follows  because 
ij  and  m are  arbitrary.  Eq.  (3.13)  is  evident  from  an  examination  of 
the  form  of  the  Green’s  function,  in  eq.  (3.12).  However,  the  symmetry 
property  is  much  more  general  than  may  be  deduced  in  this  manner, 
because  (3.12)  is  valid  only  for  free  space. 

The  magnetic  field  in  the  cavity  can  similarly  be  written  as 
2 , 2 

«(Q,k)=  / r (Q,P,k)K  (P)  dsp.  (3.15) 

2 

The  cavity  Green’s  function,  T , is  a separate,  extensive  topic,  and  we 

2 

shall  not  use  it  explicitely.  The  discussion  of  F is  therefore  deferred 
until  Appendix  II. 

* Eq.  (2.21)  has  been  derived  for  the  electric  field  of  an  electric  current  element;  obviously,  a similar 
statement  holds  for  the  magnetic  field  of  a magnetic  current. 


j 
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If  we  set  the  tangential  magnetic  fields  equal  on  the  two  sides  of 


the  aperture,  we  obtain  an  integral  equation  for  the  electric  field: 

„ f ^ (Q,P,k)K  (P)  dsp  = / r (Q.P.k)K(P)  dsp.  (3.16) 

J \ •'a 


To  get  the  invariant  formula,  multiply  by  K and  integrate  again. 


J J K (Q)  r (Q,P,k)K  (P)  dsp  dsQ 


AQ  AP 


= f f K(Q)  r (Q,P,k)K  (P)  dsp  dsQ.  (3.17) 

J A J A 


AQ  'AP 


In  Appendix  III  it  is  shown  that  values  of  k computed  from  this  formula 
are  stationary  with  respect  to  first-order  variations  of  K about  the 
true  value. 

The  stationary  property  makes  this  a good  approximation  method, 
because,  if  the  first  guess  to  the  field  is  close  to  the  true  value,  then 
the  calculated  k would  be  expected  to  be  closer  to  its  true  value,  since 
the  error  in  it  is  proportional  to  the  square  of  the  error  in  E.  This  is 
a very  rough  statement,  but  qualitatively  it  gives  the  basic  reason  why 
the  variational  method  is  sometimes  superior  to  other  approximation 
schemes . 

The  tangential  electric  field  in  the  aperture  can  be  expanded  in 
the  waveguide  modes  described  in  section  II,  and  the  expansion  coef- 
ficients and  eigenvalues  can  be  found  by  virtue  of  the  stationary 
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property  (Ritz  method).  Then,  after  these  are  known,  the  fields  every- 
where may  be  found  by  eqs.  (3.10)  and  (3.15).  Thus,  formally,  the  com- 
plete solution  of  the  problem  may  be  found.  Practically,  of  course,  the 
modes  are  known  only  for  certain  simple  shapes,  and  each  step  in  the 
Ritz  procedure  involves  very  laborious  calculations.  Hence  the  first 
step  in  this  procedure,  the  first-order  calculation  of  k,  is  as  much  as 
is  generally  done.  In  this  paper  we  go|  only  a little  further,  to  indicate 
that  one  waveguide  mode  is  a good  approximation. 

Consider  the  integral  in  eq.  (3.17): 


~ 1 

K (Q)  r (Q,P,k)  K(P)  ds 


P 


1 

H ■ z ds. 


(3.18) 


This  has  the  dimensions  of  power,  and  if  we  divide  it  by  f E^  ds  we 

JA 

have  an  admittance.  Actually,  a similar  aperture  admittance  has  been 
defined  for  a waveguide.  If  an  infinite  waveguide  radiates  through  an 
aperture  in  a ground  plane,  the  aperture  admittance,  for  the  propagating 
waveguide  mode,  is  just 


Y =- 


Il(0) 

V^O) 


(3.19) 


where  the  symbols  are  defined  in  section  II.  It  may  be  shown  that 

this  is  equal  to  3 1 

J f K(Q)[  r (Q,P)  - r (Q,P)]  K(P)  dspdsQ 


Y = 


'A  J A 

aq  ap 


f E‘  ej  ds 
-■'a 


(3.20) 


* See  any  book  on  quantum  mechanics  for  a discussion  of  this  technique. 
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where  ej  = transverse  electric  field  for  the  propagating  mode,  and 

3 00 

r (Q,P)=-  2 Yv  \>  (Q)  hv  (P).  (3.21) 

v =2 
3 

r is  written  as  a dyadic.  yv  is  the  characteristic  admittance  for  the 
vth  waveguide  mode  (yv  = l/zv,  eq.  (2.47)).  This  formula  for  Y is 
stationary  in  the  sense  that  its  first  variation  is  zero  for  variations  of 

E about  the  true  value.  It  may  be  obtained  by  equating  J E x H • z ds 

.A.  ~ 

across  a surface  consisting  of  the  aperture  and  a waveguide  cross- 
section  so  remote  that  only  the  propagating  mode  exists  there. 

It  is  also  convenient  to  define  an  “interaction  admittance”: 

YVfi  = f f hv(Q)  r (Q,P)  fa  fj.  (P)  dsp  dsQ.  (3.22) 

aq  ap 

Now  the  first  approximation  in  solving  the  problem  is  to  assume 
that  the  transverse  electric  field  in  the  aperture  is  just  the  first  mode, 
e^.  Then  we  have  immediately 

f E x H • z ds  - Y-  Yn.  (3.23) 

•'A  — — - 


For  the  infinite  waveguide  case,  the  correct  values  for  these  three 
quantities  are  actually  very  close.  Calculations  of  Yjj  have  been  made 
for  a series  of  rectangular  waveguides  for  the  steady-state  case,  and 
they  agreed  with  measured  values  of  Y to  within  experimental  ac- 

17 

curacy. 
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For  a general  cavity  we  could  still  assume  that  the  aperture  field 

was  the  first  waveguide  mode  appropriate  to  the  aperture  shape.  Values 

of  /.  Ex  H • z ds  would  then  be  available.  For  simple  geometries 
A 

Yji  may  be  calculated,  or  Y may  be  measured  for  any  shape  waveguide. 
The  measurement  would,  of  course,  be  of  a steady-state  impedance, 
but,  as  shown  below,  this  is  all  that  is  needed  for  one  approximation 
method  of  actually  performing  the  numerical  calculations. 

It  is  to  be  observed  that  the  general  procedures  used  in  setting  up 
this  problem  are  applicable  to  a much  wider  variety  of  normal-mode 
problems.  Suppose  it  is  desired  to  find  the  normal  modes  about  an 
arbitrary  finite  conducting  bofly.  The  equivalent  currents  to  be  used 
are  electric  currents  J = n * H and  K = E * n,  flowing  on  the  body  sur- 
face and  radiating  into  free  space.  These  currents  produce  the  correct 
field  exterior  to  the  body.  Note  though  that  K = 0 because  the  body  is 
a conductor.  Hence  the  currents  J in  free  space  by  themselves  pro- 
duce the  fields.  (The  currents  J are  the  actual  currents  flowing  on  the 
body.)  If  E is  the  electric  field  computed  from  J = n x H then  the 
equation 


/ 


E x H ' n ds 


0 


(3.24) 


is  true  for  the  correct  H,  and  values  of  k computed  from  this  formula 
are  stationary  with  respect  to  first-order  variations  of  H about  the 
true  value. 
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Therefore  it  is  seen  that  the  calculation  of  Yjj  for  a square 
aperture  is  exactly  the  calculation  needed  to  find  the  normal  modes  of 
a square  plate.  Also,  it  is  the  calculation  needed  to  find  the  normal 
modes  of  a square  aperture  in  an  infinite  conducting  screen,  with  free 
space  on  both  sides  of  the  screen.  These  two  systems  would  have  a 
very  low  Q. 

2.  WAVEGUIDE  CAVITY 

Now  restrict  the  cavity  to  be  a section  of  waveguide,  as  in  Fig.  13. 

The  cross  section  is  arbitrary  but  uniform 
along  the  z-axis.  The  direct  connection 
between  the  arbitrary  cavity  expansion  and 
the  waveguide  case  is  presented  in  Appendix  II; 
in  this  section  we  shall  proceed  on  the  basis 
of  waveguide  theory,  as  outlined  in  section  II. 

In  the  waveguide  region,  the  transverse 
electric  and  magnetic  fields  of  the  normal 
mode  may  be  written  as  a sum  of  waveguide  modes: 


E (xyz)  = I Vy  (z)  ev  (xy) 


H(xyz)  = X Iy  (z)  hy  (xy). 


(3.25) 

(3.26) 


Inasmuch  as  V and  I obey  transmission  line  equations,  we  have 


Iy  («fc ) = - V„  (t ) Yv  cothy  v i 


(3.27) 
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The  minus  sign  comes  from  the  fact  that  waves  traveling  in  the  positive 
z-direction  are  the  reflected  waves  from  the  short  circuit.  (See  dis- 
cussion in  section  II  above.)  Eq.  (3.27)  may  be  derived  by  combining 
eqs.  (2.35),  (2.36),  (2.44),  (2.47),  and  (2.48)  with  the  requirement  that 
111  = 0 at  z = 0. 

Thus  the  magnetic  field  in  the  aperture  becomes 

H (xyl)  = £ - yv  cothyy  l hv  (xy)  J E • ey  ds,  (3.28) 

* A 


upon  substitution  of  (3.27)  and  (3.25)  into  (3.26).  This  may  be  written  as 

t r 2 

ii(Q,k)  = j r (Q.P.k)K(P)  dsp,  (3.29) 

J A 

2 

where  K = z x E,  and  the  waveguide -cavity  Green’s  function,  r , is 


2 

r (Q,P,k)  = 2 ~ Yv  coth  yv  £ hv  (Q)  b.v  (P). 


(3.30) 


Note  that  points  P and  Q are  both  in  the  aperture. 

To  get  the  variational  expression,  multiply  (3.28)  by  z x E,  inte- 
C 1 . 

grate,  and  set  equal  to  IE  x H • zds.  This  gives 

2-y^coth  (f  E • ev  ds)2  = f f K (Q)  r( Q,P,k)  K (P)  dsp  ds 

V JAq  JAp 

(3.31) 

This  formula  is  the  invariant  expression  for  the  propagation  constant 
k.  It  is  made  more  concise  by  inserting  the  expansion  for  E,  eq.  (3.25), 
for  convenience  omitting  the  argument  z = % . This  gives 
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(3.32) 


^ ~ Yu  coth  Jv  't  Vv  = 2 2 \£,  Yv  ^ . 


As  a first  approximation,  let  all  the  V’s  be  zero  except  Vj. 
Physically,  this  means  that  the  aperture  electric  field  is  equal  to  the 
lowest  waveguide  mode  configuration.  It  is  known  that  this  cannot  be 
correct,  but  it  is  a good  approximation,  at  least  for  those  normal  modes 
which  essentially  consist  of  only  that  one  waveguide  mode.  This  ap- 
proximation leaves 


- y j coth  y j 't  = Yj  j . 


(3.33) 


This  formula  is  analogous  to  the  formula  used  in  section  I for  finding 
the  normal  modes  of  an  elementary  circuit.  If  we  consider  terminals 
at  the  aperture  for  the  lowest-order  waveguide  mode,  then 

Yleft  = ' Yright 

gives  eq.  (3.33). 

Let  us  examine  the  situation  which  results  if  two  modes  exist  at 


the  aperture.  Let 


E - Vi  e i + V2  e2. 


(3.34) 


Then  eq.  (3.32)  becomes 

2 ■ 2 „ 2 2 * 
Vj  yx  coth  7 j l + V2  y2  coth  y 2 l + Vj  Yn  + V2  Yzz  + 2V1VZ  Y1Z  = 0.  * 

(3.35) 


By  the  reciprocity  theorem,  = Ygj. 
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Call  this  last  expression  F.  Since  we  have  set  it  identically  equal  to 
zero,  it^s  stationary,  and 


and 


Now 


dF 

dV, 


= 0 


dF 

= 0 

dV2 


dF 

dV: 


B F B F dk 
3V1  +TiT  dVj  ' 


(3.36) 


(3.37) 


But  k is  the  stationary  parameter,  and  dk/dVj  = 0.  Hence 


Vj  (yj  coth  y 1 "t  + Yn)  + V2  Y12=0 

V1  Y12  + V2^2  coth  y 2 ^ + Y22)  = °- 


(3.38) 


This  is  a pair  of  simultaneous  homogeneous  equations;  the  only  way 
they  can  be  satisfied  is  for  the  determinant  to  be  zero. 


(yj  coth  y1l  + Y1;l) 
Y12 


Y 


12 


(y2  coth 7 + Yzz) 


= 0.  (3.39) 


The  expansion  of  this  determinant  provides  an  equation  for  the  eigen- 
value k . 

Consider  the  ratio  of  the  amplitudes  of  the  modes.  From  (3.38), 


^2  - Y! 2 
V1  yz  coth  7 2 1 + Y22 
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Suppose  the  first  waveguide  mode  is  “above  cutoff”  and  the  second  is 
not.  Then,  according  to  the  ideas  developed  in  section  II,  y j is  essen- 
tially imaginary  and  y 2 is  essentially  real.  The  hyperbolic  cotangent 

V2 

is  large  when  its  argument  is  small,  so  — — ♦ 0 as  ■i  - 0.  This  would 
be  expected  from  the  development  in  section  II.  Refer  to  Fig.  10,  which 
portrays  the  transverse  electric  field  of  a cutoff  mode.  As  % decreases 
there  is  less  attenuation,  and  the  total  aperture  field  decreases  in  re- 
lation to  the  “excitation.” 

The  admittances  Yj2  and  Y22  would  have  to  be  calculated  to  find 
a more  precise  effect  of  the  higher  mode.  It  is  felt  that  this  slight 
correction  is  not  worth  the  labor  involved  in  the  additional  calculation. 

All  the  calculations  made  on  the  basis  of  one  waveguide  mode  agree 
well  enough  with  experiment  to  support  this  view. 

3.  APPLICATION  TO  SQUARE  WAVEGUIDE 
We  shall  now  apply  the  first  approximation  of  the  variational 

method  to  a square  waveguide,  see  Fig.  14. 
It  is  assumed  that  only  the  TEjq  mode 
exists  in  the  waveguide.  The  electric 
field  for  this  mode  has  only  one  com- 
ponent, and  is  sinusoida-lly  distributed. 

The  desired  eigenvalues  are  the 
roots,  pn,  of  the  equation 

46 


Fig.  14.  Square  waveguide  cavity. 
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f(p) = yi  coth  y + Yll  = 0* 

(3.41) 

As  defined  above, 

p is  the  complex  frequency, 

p = co  + j £ 

(3.42) 

and 

y\ 

yi=jP^ 

(3.43) 

(3.44) 

«) 

II 

H 

(3.45) 

= p^  \±  e . 

(3.46) 

The  aperture  impedance  has  been  computed17  and  is  available  as 

a series  in  (-y^r).  When  p is  real  this  variable  is  just  (a/\),  and  for  this 
case  it  has  been  drawn  as  a graph,  in  Figs.  15  and  16. 


For  convenience,  let 


and 


aP 

2ttc 


Yu(x)  = G(x)  + j B(x). 


(3.47) 


(3.48) 


Suppose  u = ijl  o and  e = e r e D.  Then 


(3.49) 

(3.50) 

(3.51) 
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Fig.  15.  Aperture  conductance  for  a square  waveguide 


16.  Aperture  susceptance  for  a square  waveguide. 


In  terms  of  the  dimensionless  variables  x,  i /a,  and  e r,  the  eigen- 
value equation  becomes 

f(x'  i’  £ r)  = - ^J2x~  COt  [*  T V4er  x2  “ l]  + G(x)  + j B(x)  = 0. 

(3.52) 

If  real  and  imaginary  parts  of  (3.52)  are  separated,  there  results 
a pair  of  simultaneous  transcendental  equations.  Various  schemes  may 
be  used  to  find  the  solutions;  we  will  use  one  that  quickly  yields  good 
approximations  to  them.  This  problem  is  almost  identical  to  the  prob- 
lem of  finding  the  normal  modes  of  a transmission  line  terminated 
with  some  load,  as  in  Fig.  17.  In  such  a case,  when  Q > > 1,  the  voltage 
is  distributed  essentially  as  if  G were  absent.  The  major  function  of 
G is  to  contribute  to  the  loss,  not  to  the  field  distribution.  This  idea 
will  be  carried  over  to  the  cavity  problem. 


B 


T 

I 


Fig.  17.  Transmission  line. 


Let  xon  be  the  roots  of 

-i£5 


,/ 4 e r X* 
2x 


1 l / 7 

— cot  TT  — y 4 e r X - 1 +j  B(x)  = 0. 


(3.53) 


These  are  real  numbers,  xon  = (^)Qn-  Now  assume  that  the  roots,  xn, 
of  eq.  (3.52)  are  close  to  xDn,  and  expand  f(x)  in  a Taylor’s  series  about 
the  point  xQn. 
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4 


f(xn)  f(xon)  + “ xon)  + 2 ^ 2 (xn  ” xon)  + • • • ~ 

o x „ ‘ 3x‘  v 


(3.54) 


If  we  keep  only  the  linear  terms,  we  obtain,  upon  substitution  of 


G(xon)  for  f(xon)> 


xn  = xon  +' 


-G(xon) 


(3.55) 


The  derivative,  evaluated  at  xQn,  is 


^on 


(f)  l\  ,a>2  ‘ 

X on  [4  e r (Y)  Qn  - 1 


eo  l 2 
+ j 2tt  e _ / esc 

rVMn  a 


/ 72 1 3YU 

'*  «r  (f)  - 1 +-r^ 

^ on  J ox 


(3.56) 

• The,  derivative  of  is  evaluated  from  the  series  shown  in  Figs.  15 

J Q 3 B 

and  16.  The  real  and  imaginary  parts  are  just  ^ ^ and  ^ ^ y—  , 
and  are  shown  in  Fig.  18. 

The  above  equations  give  the  first  approximation  to  the  eigen- 
values xn,  and  from  these  are  obtained  Qn.  Since 


(3.57) 


then 


n 2^n  a£n 


(3.58) 
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The  real  part  of  xn  gives  the  resonance  wavelength,  in  terms  of  the 
aperture,  and  the  imaginary  part  gives  the  decay  constant  in  terms  of 
the  aperture.  Qn,  of  course,  is  a relative  measure  of  how  rapidly  the 
oscillation  dies  out. 

The  validity  of  keeping  only  the  linear  terms  was  checked  by 
computing  the  second  derivative  at  one  point.  The  point  chosen  was 
(^)  = 0.60,  e r = 1 ; this  is  a point  of  very  low  Q,  and  the  approximation 

should  be  worst  here.  The  inclusion  of  the  second  derivative  lowered 
Q by  3%  and  lowered  (y)j  by  0.7%.  This  is  only  a minor  change  and 
therefore  it  is  felt  that  the  linear  terms  alone  give  sufficiently  ac- 
curate results. 

Another  approximation  which  brings  out  a close  circuit  analogy 

B c* 

may  be  made  when  Q is  high.  From  eq.  (3.56)  it  is  seen  that 


is  the  only  term  of 


B f 


“on 


which  is  real,  and  hence  it  provides  the 


“on 


only  correction  to  the  real  part  of  xn;  i.  e.,  to  the  oscillation  frequency. 
The  calculations  for  the  first  mode  indicate  that  this  correction  is  less 
than  1%  for  Q > 7.  Hence,  for  a high-Q  cavity,  we  have  the  approxi- 
mate formulas, 


(f)  =(f) 

a.  n on 


(3.59) 


Q n=  — ^ 
n 2G  B x 


60  n B B 
2G  3 co 


(3.60) 


on 
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B is  the  total  admittance  at  the  aperture.  Eq.  (3.59)  states  that  we  may 
ignore  the  aperture  conductance;  i.  e.,  the  radiated  energy,  in  com- 
puting the  oscillation  frequency.  From  (3.60)  we  see  that  3G/3w  may 

be  ignored  in  computing  Q,  but,  of  course,  G itself  cannot  be  ignored. 

16 

Eq.  (3.60)  is  an  approximate  circuit  equation  for  Q, 

Another  problem  might  b'e  to  find  what  geometry  or  geometries 
will  resonate  at  a given  L-)  . To  be  more  general,  we  may  say  that, 

' A.  'H 

in  the  n*"*1  mode,  the  cavity  is  characterized  by  the  four  dimensionless 
l a 

parameters,—,  e r,  (y)  , and  Qn.  In  addition,  there  is  the  parameter 
a 

(y)  which  is  an  auxiliary  quantity  used  in  computing.  These  param- 
K on 

eters  are  connected  by  the  three  equations,  (3.53),  (3.55),  and  (3*58). 
Hence,  we  might  expect  that  the  arbitrary  specification  of  any  two  of 

them  would  fix  the  values  of  the  others.  Actually,  this  does  not  happen. 

I a 

Any  pair  of  values  of  (— , e r)  is  admissible  and  leads  to  unique  (y)n  and 

Qn.  This  is  obvious — the  geometry  must  specify  the  normal  modes. 

But  not  all  of  the  other  pairs  of  values  are  realizeable. 

a.  CALCULATIONS  FOB  THE  FIRST  NORMAL  MODE 

A series  of  calculations  which  brings  out  the  relations  among  the 

cavity  parameters  has  been  made  for  the  first  normal  mode.  The 

results  are  displayed  in  the  form  of  graphs  where  the  independent 

variables  are  y and  e r.  The  curves  could  have  been  drawn  in  other 

ways,  but  all  the  information  is  presented  here. 

I a 

Fig.  19  shows  the  dependence  of  — on  y for  various  values  of  e r. 
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a • 

The  curves  asymptotically  approach  that  value  of  y which  corresponds 
to  cutoff.  It  will  be  observed  that  very  high  values  of  dielectric  con- 
stant are  needed  to  resonate  the  cavity  when  y small. 

The  values  of  Q which  obtain  for  the  various  configurations  of 
Fig.  19  are  shown  in  Fig.  20.  One  conclusion  is  evident — various 
values  of  e r may  be  used  to  produce  the  same  but  one  these  values 
produces  a minimum  Q.  The  envelope  of  the  family  of  curves  in  Fig.  20 
corresponds  to  the  minimum  Q.  It  increases  very  rapidly  as  — gets 
smaller. 

The  minimum  Q may  be  separately  evaluated  by  maximizing  the 
denominator  of  Qn  (eq.  (3.58))  with  respect  to  e r and—,  and  at  the 
same  time  requiring  eq.  (3.53)  to  be  maintained.  This  may  be  done  by 
using  a Lagrange  multiplier.  Let 


A - ,5m 


-G 

3f 


17  2 + £ ^ 


(3.61) 


C = 


r 


1 


cot 


= 0, 

(3.62) 


where 


3 f 
3 x 


v + j S . 


(3.63) 


A is  the  denominator  of  Qn,  and  is  to  be  maximized.  C is  the  side 
condition  on  the  variables.  Then 
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among  cavity  parameters  for  first  normal  mode  of  square  radiating  cavity. 
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0>  00 


cavity, 


(3.64) 


of  — only.  Hence 
X 


B A 

x 3 C 
+ x _ = 

B £ r 

B e r 

BA 

B C 

Q 

+ X s - 

3(|) 

5 © 

t l2L 
B £ . 

_ B n ... 

r 

B i 

B C 

3 6 

3 6 

I 

r 

B $ 

B C 

B(i) 

a 

=4> 

A glance  at  (3.56)  reveals  that  —2-  = — = 0,  since  Yjj  is  a function 

B £ -n  A V ~ ■ - 


(3.65) 


Calculation  of  the  derivatives  leads  to  the  condition 


l 

a 


1 -^r$2 

x on. 


We  „ ^ on 


B sin^ 


*l/4er(T)  -1 


on 


4 er$  *1 


1 + 2"r  I 

X on  L 


on 


1 ~ 4ir/^T  B 

1 e Q a.  X on 


= 0. 


(3.66) 


The  simultaneous  solution  of  (3.6fe)  and  (3.62)  gives  the  values  of  e r 

l a a 

and  — which  go  with  any  (— ) ; and  with  these  (-r-)  and  the  minimum  Q 

a X on  x n 

may  be  computed  as  before,  from  eqs.  (3.55),  (3.56),  and  (3.58). 

These  calculations  were  carried  out  for  the  first  normal  mode  of 

the  square  cavity.  The  minimum  Q is  shown  in  Fig.  21  as  a function  of 
a * 

— , and  the  accompaning  parameters  — and  e r are  shown  as  functions  of 

3. 

— in  Figs.  22  and  23.  It  will  be  observed  that  the  graph  of  Q is  almost 
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Fig.  21.  Minimum  Q vs.  a/Xfor  first  normal  mode  of  square  radiating  cavity. 
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Fig.  23.  Dielectric  constant  to  obtain  minimum  (J. 
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a straight  line.  The  matching  straight  line  is  shown  dashed  in  Fig.  21, 
and  the  following  formula  is  very  accurate  for  -y  < 0.35: 

°mi„  = 0AZ4  ff)‘3  • -(3.67) 

It  is  not  known  whether  there  is  any  special  significance  to  the  cubic 
dependence  on  frequency.  At  any  rate,  the  rapid  increase  of  minimum 
Q with  decreasing  y is  in  accord  with  all  experimental  evidence. 

From  Fig.  22  it  is  seen  that  relatively  shallow  cavities  are  needed 
to  produce  the  minimum  Q,  and  that  — changes  very  slowly  with  fre- 
quency. A cavity  on  the  order  of  one-third  to  one-half  as  deep  as  it  is 
wide  will  have  a low  Q. 

Fig.  23  gives  the  value  of  dielectric  constant  to  be  used  to  obtain 
the  minimum  Q at  some  value  of  This  value  is  not  critical,  since 
Q is  a relatively  slow  function  of  e , near  the  minimum  Q.  A typical 
curve  of  Q vs.  e r for  constant  is  shown  in  Fig.  24. 
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The  dielectric  constant  for  minimum  Q increases  rapidly  with 
decreasing  -p  and  it  quickly  attains  impractical  values.  For  efficient 
operation  at,  say,  — = 0.1,  it  would  be  necessary  to  resort  to  metal 

A. 

loading  of  some  kind  (e.g.,  longitudinal  ridges,  etc.).  Fundamentally, 
then,  the  cavity  would  have  a different  shape;  it  would  no  longer  be  a 
square  waveguide  cavity,  dielectric  loaded.  The  above  curves  and  dis- 
cussion of  minimum  Q were  derived  for  a square  waveguide  cavity,  but 
it  is  felt  that  they  give  at  least  the  order  of  magnitude  of  what  could  be 
expected  from  other  configurations.  On  this  basis  then  we  would  say 
that  the  minimum  Q obtainable  from  a dielectric  - loaded  square  wave- 
guide cavity  with  ~ = 0.1,  is  Q = 424.  If  by  necessity  or  choice  we  went 
to  some  other  configuration,  with  the  same  y,  we  would  not  expect  to  ap- 
preciably lower  the  Q. 

b.  EFFECT  OF  DIELECTRIC  LOSS 
The  effect  of  a finite  conductivity  of  the  medium  in  the  cavity  may 
be  found  by  using  a complex  dielectric  constant.  To  see  this,  go  back 
to  one  of  Maxwell’s  curl  equations. 


V x 


where 


Pin  = j Pn  e En  + o'  En  = J Pn  e'  En, 


e'  = e (1  + ) 

JPn6' 


(3.68) 


- 6 (*  + — 7 ) = e (1  - j tan  S ).  (3.69) 


m e 

n 


This  approximation  is  valid  when  Qn  > > 1.  The  complex  relative 


■srr«iir.g<iw-r- 
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dielectric  constant  may  be  introduced: 


where 


e*  = e o e;  * (3-7°) 

e;  = er  (1  - j tans  ).  (3.71) 


This  complex  relative  dielectric  constant,  e may  be  inserted  into  the 
above  equations  in  place  of  ^r,  and  the  real  and  imaginary  parts  re- 
evaluated. Rather  than  do  this,  it  is  easier,  since  tan  S < < 1 for  prac- 
tical dielectrics,  to  consider  f (eq.  (3.52))  as  a function  of  two  complex 
variables.  Then  f may  be  expanded  in  a Taylor’s  series,  and  only  the 
linear  terms  need  be  retained.  Thus, 


f(xn,  e ’ ) = f(xon,  g r)  + 

3 x 


(xn  - xon)  + 


xon  er 


3 f 
Be 


on. 


The  second  term. 


3 f 
3 x 


on. 


evaluated  from  eq.  (3.52). 


(3.72) 

, is  given  in  eq.  (3.56).  The  third  term  is 
e 

r 

The  result  is  (after  substitution  of  (3.53)), 


G = 


3 f 
3 e ' 


(e  r “ e r) 


xon*  e r 


= •]  rrj—^  — (— ) esc2  it Ja  e 

1 ^ a ^ on  L a V * 


(r>2  - 1 

K on 


2(f)__  B 


on 


4e  r(r")  -l 

r K on 


From  eq.  (3.72),  the  eigenvalues  are  given  by 

G + 0 


*ns*on-|7 

3 x 


-e  r tan  § 


(3.73) 


(3.74) 


xon»e  r 
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G is  the  aperture  conductance;  we  may  call  G the  radiation  conductance 


and  (J  the  dielectric  conductance.  When  the  Q is  high, 
imaginary,  and 


I xon>  er 


■ <T), 


(3.75) 


^n 


3 x xon»  6 r 
2(G  + Q.  ) 


(3.76) 


The  dielectric  loss  does  not  affect  the  resonance  frequency  so  long  as 
Q remains  high,  but  it  may  have  an  appreciable  effect  on  Q. 

From  eq.  (3.76)  one  may  write 


1 11 

Q™  Q Q ' 
T D 


(3.77) 


where  Qjq  is  the  radiation  Q,  Qj~,  is  the  dielectric  Q,  and  Qq-  is  the 


total  Q.  Then 


Qn  3 1 

D ^Tjn  3 x 


(3.78) 


Now  numerical  calculations  for  the  first  mode  reveal  that  the  csc^  terms 

in  Q and  are  at  least  five  times  as  big  as  the  others.  If  they  are 

dropped,  -i—  immediately  reduces  to  tan  S . Thus 
QD 


1 1 


+ tan  S . 


(3.79) 


A “good”  dielectric  has  tan  S on  the  order  of  0.01  or  smaller.  This 
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gives  an  indication  of  the  effect  of  the  dielectric  on  the  total  Q. 

A few  statements  about  the  stored  energy  associated  with  the 
radiating  cavity  may  be  made  by  referring  to  eq.  (1.12),  the  energy 
definition  of  Q.  By  comparing  this  definition  with  eq.  (3.76)  it  is  seen 
that 


3 f 


3 x 


must  be  a measure  of  the  stored  energy;  G,  a measure 

|xon»  e r 

of  the  radiated  energy;  and  Q,  a measure  of  the  energy  dissipated  in 


the  dielectric. 

Refer  now  to  eq.  (3.56).  The  individual  terms  may  be  ascribed  to 
contributions  to  stored  energy  from  different  regions  of  the  total 
volume.  Thus,  the  first  two  terms  are  interpreted  as  a measure  of 


the  energy  in  the  waveguide,  and  ^ 

<3 


3 B 


is  a measure  of  the  energy 


'-on 


in  the  half-space.  Now,  as  noted  above,  the  csc^  term  is  much  bigger 
than  the  others,  so  we  may  say  that  most  of  the  stored  energy  is  in  the 
cavity,  and  very  little  is  in  the  half-space. 

We  can  now  find  a more  physical  reason  why  the  reciprocal  of  the 
dielectric  Q equals  tan  8 . A closed  cavity  with  lossless  walls  has 
= tan  S . If  the  radiating  cavity  has  a high  Q most  of  the  stored 
energy  will  be  inside  the  cavity.  Also,  the  fields  in  the  cavity  will  be 
close  to  the  configuration  existing  in  a closed  cavity.  Since  the  stored 
energy  and  losses  approximate  those  in  a closed  cavity,  it  is  expected 
that  tv—  = fan  § . 

* For  an  interesting  discussion  on  location  of  electromagnetic  energy,  see  Mason  and  Weaver, 

The  Electromagnetic  Field,  Dover  Publishing  Company,  New  York,  1929.  p.  266. 
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The  radiation  efficiency  may  be  found  from  these  energy  con- 
siderations. From  the  statements  concerning  G and  Q , it  is  seen  that 


Efficiency  = ^ , (3.80) 


which  may  also  be  written  as 

Q-p 

Efficiency  = . (3  81) 

Qr 

If  a substitution  is  made  from  eq.  (3.79),  one  has 


Efficiency  = . (3.82) 

1 + Qr  tan  8 


This  formula  illustrates  the  principle  that  a high-Q  system  fundamentally 
has  a low  efficiency.  Physically,  this  is  due  to  the  high  stored  energy; 
i.  e.,  to  high  circulating  currents  and  high  local  field  intensities. 
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IV.  CORRELATION  OF  NORMAL  MODES  WITH 
STEADY-STATE  BEHAVIOR 

1.  THEORY 

The  correlation  of  normal  modes  of  a linear  system  with  steady- 
state  behavior  is  a subject  which  has  been  widely  discussed,  particu- 

19 

larly  in  reference  to  circuits.  Foster  has  shown  that  any  finite  loss- 
less network  is  equivalent  to  a chain  of  elementary  LC  circuits  each 
of  which  is  resonant  at  one  of  the  normal  fnode  frequencies,  and  has 
stored  energy  corresponding  to  that  in  the  normal  mode.  The  theorem 
is  readily  extended  to  cover  slightly  lossy  networks  by  the'addition  of 

an  R into  each  LC  circuit.  This  is  an  approximation  which  is  valid 

2 0 

for  high-Q  networks.  Schelkunoff  has  discussed  this,  and  the  exten- 
sion to  a system  with  an  infinite  number  of  modes. 

The  cavity  antennas  considered  in  this  paper  are  in  the  class  of 
slightly  lossy  networks  with  an  infinite  number  of  degrees  of  freedom. 
Hence,  it  is  assumed  that  the  antenna  is  representable  by  a chain  of 
RLC  circuits,  as  in  Fig.  25.  The  antenna  is  connected  to  a transmis- 
sion line  of  characteristic  impedance  ZQ. 

The  application  of  this  equivalent  circuit  (and  the  alternate  one, 
a chain  of  series  circuits)  to  cavity  problems  has  been  discussed  in 

the  M.I.T.  Radiation  Laboratory  Series,21’22  and  by  Slater.23’24  Some 

• * 

of  the  results  we  shall  use  will  be  in  a slightly  different  form  from 
those  in  the  references,  but  the  derivations  are  similar  and  will  not  be 
repeated  here. 
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Any  one  of  the  elementary  RLC  circuits  of  Fig.  25  has  an  imped- 


i* 


ance  which  changes  rapidly  near  resonance,  but  is  very  stable  away 
from  resonance.  If  one  of  these  circuits  is  resonant  at  a frequency  far 
from  that  of  any  other  circuit,  then,  in  the  vicinity  of  its  resonance,  the 
rest  of  the  network  may  be  lumped  into  the  series  impedance  Za,  which 
is  essentially  constant.  The  resulting  equivalent  circuit  is  shown  in 
Fig.  26.  It  is  now  assumed  that  the  normal-mode  Q is  related  to  the 
circuit  by 


Qn  = 


wnL 


(4.1) 


and  that  the  normal-mode  frequency  is  equal  to  the  resonance  frequency, 

^ = -— = . In  the  references  cited  this  Q is  spoken  of  as  the  unloaded 
VLC 

Q,  unloaded  in  the  sense  that  ZQ  is  not  included. 

A curve  of  standing  wave  ratio  (SWR)  versus  frequency,  for  this 
circuit,  would  have  a resonance  dip,  and  the  sharpness  of  it  would  de- 
pend on  Q.  Typical  curves  are  shown  later  (see,  for  example,  Fig.  29). 
This  is  a measured  curve  for  an  antenna,  but  it  is  assumed  that  there 
is  an  elementary  circuit  which  matches  it.  The  frequency  of  minimum 
SWR  is  the  resonance  frequency. 

The  bandwidth  is  a term  which  is  often  used  to  describe  the  oper- 
ating range  of  an  antenna  in  terms  of  some  specified  SWR  limit.  How- 
ever, for  simplicity,  we  shall  define  bandwidth  as  the  reciprocal  of  the 
Q,  and  the  SWR  limit  will  be  taken  as  the  variable. 


i * 


i \ 


% 

% 
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There  are  two  cases  to  be  considered  for  the  equivalent  circuit. 
The  simpler  one  assumes  that  the  effect  of  Za  is  insignificant  and  that 
it  can  be  ignored.  In  this  case  the  following  result  can  be  shown*  Let 


^ R 


R 


cr 


f SWR  at  resonance,  if  — > 1 
= J Z, 


1 


SWR 


at  resonance,  if  £_  < 1. 


(4.2) 


Then 


Qn  = 


"n 


"2-^1 


(4.3) 


where  ^ and  ^ are  the  frequencies  at  which  the  SWR  has  risen  to 
the  value  r j : 


rl  = 


1 +J*oZ  + ^ +4(T02  ~ a Q 


(4.4) 


a/1  + ? o-Q  + °o  ~ 'J1  + 2 aO 


Co 


Experimentally,  the  correct  value  for  crQ  is  obtained  by  examining  the 
impedance  curve,  to  determine  whether  R/Zq  is  greater  or  less  than 
one  at  resonance.  See  the  references  for  a discussion  of  this  point. 

Slater  has  examined  the  circuit  when  Za  is  included.  He  assumes 
that  Za  is  constant  over  a wide  frequency  range  and  that  the  SWR  far 
from  resonance  is  what  would  be  obtained  at  the  resonance  frequency, 
if  the  resonant  circuit  were  removed.  Under  these  conditions  it  is 


See  Slater’s  article.  He  finds  SWR  limits  for  the  loaded  Q.,  but  the  derivations  are  similar. 
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Fig.  25.  Equivalent  circuit. 


Fig.  26.  Equivalent  circuit  near  resonance. 


i 

! 

t 
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possible  to  show  that 


vA1  + °o)2  + (!  + ^l)2  +y (1  - a0)2  + (1  - ai)Z 

r2  " ~ ’ 

aA1  + o-0)2  + (1  + ^i)2  “V^1  " ^o)2  + (x  " c^)2 

(4.5) 

where  r ^ is  the  standing  wave  ratio  at  ^ and  o-  is  l/SWR  far 

from  resonance,  and  cr  is  defined  above. 

The  two  formulas,  eqs.  (4.4)  and  (4.5),  give  similar  results  when 
Q»  1.  With  their  aid  the  Q of  an  antenna  may  be  computed  from  an 
experimental  measurement  of  SWR  vs.  frequency  (plus  some  phase 
measurement  to  eliminate  the  ambiguity  about  rr0).  Conversely,  if  one 
has  a value  for  Q,  the  SWR  limits  on  the  bandwidth  may  be  computed, 
after  some  minimum  SWR  has  been  assumed.  In  the  antenna  problem 
it  may  be  assumed  that  the  minimum  SWR  is  one,  since  the  feed  geome- 
try is  arbitrary  and  can  be  changed  until  this  value  is  obtained.  In  this 
case  the  bandwidth  is  the  percentage  frequency  range  within  which  the 
SWR  is  less  than  2.62.  This  value  is  obtained  from  eq.  (4.4)  when 
^Q=l,  and  from  eq.  (4.5)  when  tq  = 1 and  cq  = 0. 

It  is  possible  to  get  lower  SWR  limits  on  the  bandwidth  than  the 
value  given  above;  namely,  SWR  = 2.62  for  minimum  SWR  =1.  In 
eq.  (4.4),  r^  may  be  minimized  with  respect  to  c rQ;  the  result  is 
r j = 2.42,  for  CTo  = -fl.  Hence,  if  crQ  =/" 2,  one  has  the  minimum  SW& 
at  the  ends  of  the  frequency  band.  This  remark  is  of  some  interest 
because  the  practical  problem  is  to  minimize  the  SWR  as  much  as  possible. 

AF  18(600)85  486-7  72 


0 


2.  EXPERIMENTAL 


Three  experimental  models  were  constructed  and  tested  to 
check  the  theory.  These  were  all  square  waveguide  cavities,  filled 
with  a dielectric,  and  flush  mounted  in  a ground  plane.  These  antennas 
were  excited  in  various  ways  and  the  input  impedance  was  measured. 
From  curves  of  SWR  vs.  frequency  the  resonance  frequency  and  Q 
were  obtained,  as  described  above. 

The  first  cavity  was  3 by  3 inches  in  cross  section  and  23  inches 
deep.  It  was  filled  with  a block  of  cast  sulfur  and  mounted  in  a 3 -foot 
square  ground  plane.  The  dielectric  constant  of  the  sulfur  was  meas- 
ured by  closing  the  aperture  with  an  aluminum  plate  and  finding  the 
resonance  frequencies  of  the  closed  cavity.  A simple  calculation  gave 
the  resonance  frequencies  of  the  air-filled  closed  cavity  (as  a check, 
these  were  also  measured),  and  the  two  sets  of  frequencies  were  con- 
nected by  ^/i^.  This  procedure  yielded  a value  for  the  dielectric  con- 
stant, £r  = 3.80,  with  an  average  deviation  of  0.01  (for  the  first  six 
resonances). 

The  antenna  was  driven  by  a loop  in  the  rear  wall,  as  shown  in 
Fig.  27.  The  short-circuited  stub  in  series  with  the  loop  was  used  to 
help  the  impedance  match.  The  exact  configuration  of  the  loop  and 
stub  is  unimportant,  since  we  are  interested  only  in  the  resonance 
frequency  and  Q.  These  characteristics  are  independent  of  the 
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excitation,  since  it  is  assumed  that,  near  resonance,  the  fields  are  es- 


sentially the  normal-mode  fields.  The  feed  geometry  controls  the 
magnitude  of  the  impedance  at  resonance,  and,  of  course,  it  does  per- 
turb the  field  and  change  the  resonance  properties  somewhat.  This 
effect  will  be  noted  later. 

The  antenna  was  mounted  on  the  side  of  the  building  and  allowed 
to  radiate  into  space.  The  input  impedance  was  measured  with  a pre- 
cision slotted  line.  It  is  estimated  that  the  measured  frequency  is 
accurate  to  within  2 megacycles,  and  that  small  frequency  intervals  (used 
in  measuring  bandwidth)  are  accurate  to  within  0.2  megacycle.  The 
values  of  SWR  are  probably  correct  to  within  5%  for  SWR  < 50.  Higher 
values  of  SWR  are  not  accurate  because  of  losses  in  the  slotted  line. 

The  SWR  of  the  antenna  is  shown  as  a function  of  frequency  in 
Fig.  28. * The  region  from  1000  to  1026  megacycles  is  shown  dashed 
because  the  theory  predicted  a resonance  dip  there.  It  was  missed  in 
the  measurement.  The  structure  of  the  other  resonance  dips  was  care- 
fully measured,  and  a typical  one  is  shown  in  Fig.  29.  The  crosses 
indicate  the  bandwidth  points;  the  outer  ones  are  at  values  of  and 
the  inner  ones  are  at  r£.  The  phase  of  the  impedance  was  also  meas- 
ured, although  it  is  not  shown.  It  turned  out  that  R/Zq  was  less  than 

* For  convenience,  all  standing  wave  ratio  curves  are  plotted  on  a logarithmic  scale. 
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Fig.  28.  Standing  wave  ratio  vs.  frequency,  for  deep  suifur-filled  antenna. 


r 


one,  so  cr Q = l/SWR  for  each  mode.  This  fact  indicates  that  there  was 
weak  coupling  between  the  transmission  line  and  the  cavity.  A better 

i 

impedance  match  could  have  been  obtained  by  using  a bigger  feed  loop. 

The  antenna  has  £/ a = 7.67,  and  the  first  normal  mode  is  very 
close  to  cutoff,  which  occurs  at  1009  megacycles.  (This  frequency  is 
indicated  in  Fig.  28.)  The  higher  modes  come  in  quick  succession  I 

i 

because  the  guide  wavelength  changes  very  rapidly  with  frequency.  The 
first  six  normal-mode  frequencies,  and  the  Q's,  were  computed  accord- 
ing to  the  above  theory.  The  loss  factor  (tan  8 ) was  taken  to  be  the  value 
given  by  Von  Hippel,25  namely,  tan  8 = 0.0025.  This  changes  slightly 
with  frequency. 

All  of  the  results  are  in  Table  I.  The  second  and  third  columns 
give  the  measured  and  calculated  resonance  frequencies  for  the  various 
modes.  The  fourth  column  gives  the  measured  Q assuming  only  a 

j 

simple  parallel  circuit,  by  eq.  (4.4).  Column  5 gives  the  measured  Q 
using  the  more  refined  circuit,  by  eq.  (4.5).  Column  6 gives  the  Q cal- 
culated from  radiation  losses  only,  while  column  7 gives  the  calculated 
total  Q.  Column  8 gives  the  efficiency. 

It  will  be  noted  that  the  calculated  and  measured  resonance  fre- 
quencies agree  very  closely.  In  most  cases  the  measured  values  are 
a trifle  higher  than  the  calculated  ones.  The  Q values,  however,  are 
not  so  close.  In  the  first  place,  there  is  a considerable  difference  be- 
tween the  radiation  Q and  the  total  Q.  This,  of  course,  is  due  to  the  fact 

i 

i 

1 
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that  the  cavity  is  very  deep  and  the  dielectric  losses  are  large.  This 
is  particularly  true  for  the  first  modes,  near  cutoff.  Very  close  to 
cutoff  there  are  high  standing  waves  in  the  guide,  and  these  cause  the 
excessive  dielectric  losses  which  result  in  a low  Q and  a low  efficiency. 

TABLE  I 

Results  for  3 by  3 by  23 -inch  sulfur-filled  loop-fed  cavity  antenna. 


Mode 

Resonance 

Frequency 

Measured  Q 

Calculated  Q 

Calculate! 

Measured 

Calculated 

Qi 

°R 

Qq1 

Efficiency 

1 

- 

1016 

- 

- 

1470 

344 

23% 

2 

1036 

1037 

104 

134 

430 

220 

51% 

3 

1074 

1072 

89.6 

114 

240 

156 

65% 

4 

1122 

1120 

82.5 

97.7 

172 

124 

72% 

5 

1184 

1180 

76.9 

90.4 

137 

105 

77% 

6 

1254 

1251 

67.1 

75.6 

119 

93.4 

79% 

The  calculated  Q^.  is  to  be  compared  with  the  measured  Q.  The 
Q2  values  are  higher  than  the  ones  and  so  are  closer  to  the  com- 
puted values;  they  would  be  expected  to  be  closer  because  they  are 
derived  from  a better  equivalent  circuit.  In  every  case  the  computed 
Q-j-  is  higher  (in  some  cases  much  higher)  than  Q£.  This  also  is  to  be 
expected,  since  metal  losses  were  ignored  in  the  calculation.  The  dis- 
crepancy is  greatest  where  Q is  very  high,  because,  in  general,  as 
noted  above,  a high-Q  system  has  high  losses  and  a low  efficiency. 
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The  second  experimental  antenna  was  a 4 by  4 by  8rinch  square 
cavity,  also  filled  with  sulfur.  The  dielectric  constant  of  the  sulfur 
was  er  = 3.88.  The  difference  in  er  for  the  two  sulfur  blocks  is  appar- 
ently due  to  the  manner  in  which  it  was  prepared  - the  cooling  rate  of 
liquid  sulfur  has  a considerable  effect  on  crystal  structure  and  air- 
bubble  formation. 

This  antenna  was  excited  by  a stub  close  to  the  aperture,  as  shown 

“ in  Fig.  30.  This  location  for  the 
stub  was  chosen  because  it  is  near 
a voltage  maximum  for  the  first 
two  normal  modes.  Fig.  31  shows 
the  magnitude  of  the  electric  field 
for  the  TE1q  waveguide  mode  for 
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Fi".  30.  Second  sulfur-filled  antenna. 


the  first  two  normal  modes.  These 


curves  were  drawn  by  finding  the  guide  wavelength  at  the  two  resonance 


of  second  sulfur-filled  antenna. 
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A good  impedance  match  to  the  slotted  line  was  obtained  with  a 
stub  1 inch  long  and  l/8-inch  in  diameter.  This  stub  screwed  on  to  a 
modified  UG-58U  fitting,  so  that  the  total  protrusion  into  the  cavity  was 
about  3/l6-inch  longer  than  the  probe.  Other  length  probes  were  also 
tried,  and  they  produced  different  resistance  values  at  resonance. 

Fig.  32  shows  the  SWR  vs  frequency  for  the  1-inch  probe,  and  Fig.  33 
shows  the  SWR  for  a 0.3-inch  probe.  It  should  be  remembered  that 
the  data  for  the  0.3-inch  probe  are  not  as  accurate  as  those  for  the 
1-inch  one,  because  of  the  high  SWR’s  involved. 

There  are  four  major  resonance  dips  present  on  the  curves  in 
Figs.  32  and  33.  The  first  two  are  due  to  normal  modes  involving  only 
the  TE^q  waveguide  mode.  These  two  are  compared  with  theory  in 
Table  II.  The  TE^j  mode  propagates  after  the  second  dip  (cutoff  fre- 
quency of  TEjj  mode  is  shown  on  graphs),  and  the  third  dip  is  due  to 
the  first  resonance  of  the  TE^  mode.  This  portion  of  the  curve  is 
like  that  obtained  from  two  coupled  circuits,  and  that  is  exactly  the 
situation  here,  since  there  are  two  propagating  modes  coupled  at  the 
probe.  The  third  resonance  of  the  TE^q  mode  is  calculated  to  be  at 
1210  megacycles,  so  this  must  correspond  to  the  fourth  dip  of  the 
curves.  Again,  it  is  a multiple  dip  because  there  are  two  propagating 
modes.  If  the  feed  were  balanced,  so  that  the  TEj^  mode  were  not 
excited,  then  the  third  dip  would  be  absent,  and  the  fourth  would  be  clean. 
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Frequency  — Megacycles 


TABLE  II. 


Results  for  first  two  normal  modes  of  4 by  4 by  8-inch 
sulfur -filled,  probe -fed  antenna. 


Probe 

Mode 

Resonance  Frequency 

Measured  Q 

Calculated  Q 

Calculate 

Measured 

Calculated 

Q1 

qr 

Qq1 

Efficient 

1 -inch 

1 

782 

799 

39.6 

41.0 

71.9 

61.9 

86% 

2 

936 

961 

30.7 

32.4 

32.7 

30.4 

93% 

3 -inch 

1 

794 

799 

35.2 

41.8 

71.9 

61.9 

86% 

2 

957 

961 

22.3 

26.0 

32.7 

30.4 

93% 

The  first  two  resonance  dips  were  measured  carefully,  and 
curves  similar  to  that  in  Fig.  29  were  obtained.  It  will  be  noted  from 
Table  II  that  the  resonance  frequencies  of  the  small  probe  are  closer 
to  the  computed  values  than  those  of  the  larger  probe.  The  small  shifts 
in  resonance  frequency  are  due  to  a perturbation  of  the  fields  by  the 
probe;  the  smaller  probe  causes  less  perturbation  and  hence  less  shift. 
According  to  general  perturbation  theory  of  cavities26  a metal  protru- 
sion in  a region  of  high  E and  low  H should  lower  the  resonance  fre- 
quency, whereas  one  is  a region  of  low  E and  high  H should  raise  it. 
This  general  statement  correctly  predicts  the  lowering  of  the  resonance 
frequency  in  this  case.  Also,  according  to  this  theory,  a loop  in  the 
rear  wall  should  raise  the  resonance  frequency,  and  this  was  observed 
in  the  first  antenna. 
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Observe  that  the  measured  Q2  for  the  1-inch  probe  at  the  second 
mode  is  higher  than  the  calculated  Qp.  This  is  the  only  case  in  all 
the  measurements  where  this  occurred,  and  it  is  not  expected  because 
wall  losses  should  lower  the  measured  Q from  the  computed  value. 
However,  this  measurement  is  unique  in  one  other  respect  — R/Zq  is 
greater  than  one,  so  that  the  transmission  line  is  strongly  coupled  to 
the  cavity.  This  strong  coupling  implies  that  there  will  be  high  local 
fields  near  the  probe,  and  apparently  the  total  stored  energy  increased 
by  more  than  enough  to  offset  the  wall  losses. 

It  was  attempted  to  make  an  independent  check  on  the  theory  by 
measuring  the  efficiency  of  the  antenna.  This  was  done  on  the  low- 
frequency  efficiency  meter,  a special  instrument  built  at  this  laboratory, 
which  is  described  in  various  reports.27  The  measured  efficiency  for 
the  1-inch  probe  at  the  first  resonance  was  about  80%.  This  is  quite 
close  to  the  computed  value,  but  it  has  tentatively  been  shown  that  the 
measured  efficiency  varies  greatly  for  different  probe  positions  and 
for  slightly  different  frequencies.  This  indicates  that  the  efficiency  is 
quite  dependent  upon  the  feed  perturbation;  i.e.,  upon  higher  nonpropagat- 
ing waveguide  modes.  This  has  only  a minor  effect  upon  the  resonance 
frequency,  but  the  effect  on  Q may  be  considerable.  Actually,  no  defi- 
nite conclusions  may  be  drawn  because  of  possible  large  inaccuracies 
in  the  efficiency  measurement. 
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The  third  antenna  was  a shallow  one  having  approximately  the 


geometry  necessary  for  minimum  Q.  It  was  basically  3 by  3 by 
1.2  inches,  but  to  avoid  difficulty  with  higher  modes  it  was  cut  in 
half  with  an  image  plane.17  This  automatically  balanced  the  antenna, 
so  that  the  TE^  mode  could  not  be  .excited.  A sketch  of  the  antenna 
is  shown  in  Fig.  34.  The  cavity  was  filled  with  a mixture  of  lead- 
chloride  and  mineral  oil,  which  had  a dielectric  constant  of  8.41  t 0.06, 
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Frequency  — Megocycles 

Fig.  35.  Standing  wave  ratio  for  shallow  half-cavity,  filled  with  lead-chloride,  mineral-oil  mixture. 


The  cavity  was  excited  with  a probe  on  the  center  line  of  the 
image  plane  and  l/4-inch  from  the  aperture.  A probe  0.240-inch 
long  and  l/8-inch  in  diameter  gave  a reasonable  impedance  match. 

The  SWE  on  the  slotted  line  was  measured  as  above  and  it  is  shown 
in  Fig.  35  as  a function  of  frequency.  The  antenna  impedance  was  also 
measured  with  a much  smaller  probe,  l/8-inch  long  and  0.05-inch  in 
diameter.  The  standing  wave  ratio  with  this  probe  is  also  shown  in 
Fig.  35. 

The  calculated  resonance  frequency  for  this  antenna  is  1200  mega- 
cycles, so  there  is  a considerable  deviation  here.  According  to  the 
above  general  theory  the  shift  is  in  the  correct  direction  and  is  less 
for  the  smaller  probe,  as  expected.  Apparently  for  shallow  cavities 
the  effect  of  the  feed  becomes  more  important.  In  this  case  it  pre- 
sumably shifted  the  resonance  frequency  by  8%.  . 

Table  III  gives  the  comparison  between  measured  and  calculated 
resonance  frequency  and  Q.  The  losses  in  the  dielectric  were  not 
considered,  because  the  value  of  tano  for  this  medium  is  not  known. 

TABLE  III. 


Results  for  first  normal  mode  of  3 by  1.5  by  1.2-inch  cavity 
with  image  plane,  filled  with  lead-chloride,  mineral  oil  mixture. 


Probe  Resonance  Frequency  Measured  Q Calculated  Q 


Large 

Small 


Measured 

Calculated 

Q! 

°2 

1100 

1200 

13.7 

15.0 

11.8  14.4 
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APPENDIX  I.  FIELDS  OF  AN  IMPULSE  CURRENT 


We  now  investigate  the  fields  produced  by  the  impulse  current 
JQ  = q 5 (P  - Q)  § (t  - tQ).  The  only  unusual  thing  about  this  current  is 
the  S -function  in  time.  The  idea  of  a current  element  which  is  in- 
finitesimal in  its  spatial  extensions  has  long  been  a convenient  concept. 
JQ  may  consist  simply  of  a point  charge  Q moved  from  A to  B at 
time  tQ.  Points  A and  B are  very  close  together,  and  the  motion 
takes  place  in  a very  small  time.  The  current  will  therefore  consist 
of  a pulse  at  t = tQ. 

There  will  be  static  fields  before  and  after  t0,  but  they  will 

differ  by  a static  dipole  field  due  to  a dipole  Q AB.  To  see  this,  use 

the  law  of  superposition.  Fig.  36  (a)  shows  the  configuration  for  t > tQ 

and  Fig.  36(b)  shows  the  initial  configuration,  for  t < tQ.  The  difference 

— ► 

between  these  is  the  dipole  Q AB.  It  is  shown  below  that  the  dipole 
moment  Q AB  must  equal  q. 


Q Q 

-Q 

Q 

A 

B A 

B A 

B 

(a)  (b)  (c) 


Fig.  36. 

Mathematically,  the  pulse  has  been  represented  by  a S -function. 
It  is  sufficient  for  our  purposes  to  define  this  function  by  its  integral, 
as  was  done  in  section  II.  For  computing  the  fields  we  will  need  the 
derivative  of  the  current,  and  thus  we  need  the  derivative  of  the 
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8 -function.  Call  this  the  D-function,  and  define  it  by 
A 


f D(r  - tQ)  dr  = S (t  - tD). 
J -00 


(A.l) 


If  the  8 -function  is  pictorially  represented  by  a steep  pulse,  then  the 


,1  D(t-t„) 


D-function  must  be  a double  pulse,  as 

in  Fig.  37.  The  D-function  is  some- 

t times  called  the  “unit  doublet  impulse" 
2 0 

function. 

Now  consider  the  equations  that 
t the  fields  satisfy.  Maxwell’s  equations 


Fig.  37.  8 and  D functions. 


3D 

Vx  H =— =+  JQ, 

O t ““ 


„ „ 3B 

7xE=-rt  • 


V . B = 0 

Ta.2) 

V • E = p . 


then 


B = V x A, 

(A.3) 

3 A 

(A. 4) 

E = - V $ 

- d t 

3$ 

(A. 5) 

A + Ate  —7  = °, 

d t 

20 

then  both  A and  <£•  satisfy  the  wave  equation,  and  A is  given  by 

A(Q)  = M / W dv  = ^ Mt  - to  - i).  (A.6) 
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By  [ Jq]  is  meant  the  retarded  value,  J0(t  - — ),  where  r is  the  dis- 
— — c 

tance  from  P to  Q.  These  equations  are  valid  for  an  isolated  source; 

_ i.  e.,  one  in  free  space. 

Now  let  q be  directed  along  the  z-axis,  at  the  origin  of  coordinates.  j 

The  coordinate  system  is  shown  in  Fig.  38.  The  vector  potential  A j 

has  only  a z-component  and,  by  symmetry,  has  no  ^ -dependence; 

therefore  H has  only  a ^-component. 

Z 


*v  =- 


_1_  3 Az  _ _ ^ Az  3 r 

pL  B p p.  B r Bp' 


(A.7) 


B r P 

Now  - — =—  = sin  0 (the  coordinate  z is  kept  constant  in  this  differ- 
3 P r 

entiation),  and,  from  (A. 6), 
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s 


— = &.R  [ _ 1 s (t  - tQ  - -)  - — D(t  - tQ  - — ) 
4it  [ . r1 2 *  ° c'  rc  v ° c 


Hence 


_ q sin  6 
0 4u 


rc  c TZ  c 


Now  find  the  potential  $ , from  eq.  (A. 5). 


B „ . 3 Az 

a e = - V • A 

B t ~d  z 


BA„  B 


■z  0 r 


3 r 3 z 


cos  0 


Substitute  (A. 8) 


B$  yu.  q cos  0 f 1 


Integrate  with  respect  to  t: 


e $ = 


q cos  0 


(A.8) 


(A. 9) 


(A. 10) 


(A.  1 1 ) 


— S (t  - tQ  -~)  U(t  - tQ  -|)  . (A.12) 

rc  c r4'  ^ 


The  integral  of  the  S -function  is  U,  the  step  function,  defined  by 


U(t  - t0)  = 


0 if  t < tr 


(A. 13) 


1 if  t > t^ 


Note  that  an  arbitrary  function  which  is  independent  of  t may  be  added 


to  3>  . This  corresponds  to  a residual  static  field. 


The  electric  field  is  now  found  from  eq.  (A. 4).  The  contribution 


from  the  vector  potential  is 


14  = 2 4^  D(t  ’ *o  - f)  = D(t  - tQ  - (r  cos  0 - 0 sin  0).  (A.14) 
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z,  _r,  6,  etc.  are  unit  vectors  along  the  respective  coordinate  lines. 

The  other  part  of  E is 

„ , B 0 1 B * 

B r — r B 0 


= r 


q cos  8 
4tr  e 


L D{t  - tQ  - (t  - tQ  - ~)  - — U(t  - tQ  - |-) 

r J 


L rc 


r c 


- 0 


q sin  8 
4tt  e 


L r c 


i_  8 (*  - t0  - Jl)  + Jj  U(t-t0--|) 


(A. 15) 


Then 


E = 2q  cos  9 /F 

r 4n 


4 8(t  - to  -7)  +^f  U(t  - t0  -i) 


l r 


c'  3 
r 


(A. 16) 


Ee  = 


q sin  8 
4ir 


rc 


D(t  - tQ  - ~)  + (t  - t„  - J)  + — U(t  - tc 


(A.  17) 

Equations  (A. 9),  (A.16),  and  (A.  17)  give  the  electric  and  magnetic 
fields  of  the  impulse  current.  In  section  II  these  fields  are  called  GQ 
and  Fq. 

The  terms  of  E and  H form  the  customary  “far  field.”  When 
the  other  terms  are  ignorable  the  field  becomes  transverse,  and,  as  in 
a plane  wave, 


(A. 18) 


The  — terms  contain  the  D-function,  so  the  field  at  a distant  point 
r 

consists  of  two  pulses  of  opposite  sign,  one  after  the  other.  They  arise 
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from  the  two  accelerations  of  the  charge  — one  to  start  it  at  A and  the 
other  to  stop  it  at  B.  This  is  analogous  to  the  close  waves  of  compres- 
sion and  rarefaction  which  propagate  after  an  explosion. 

The  -i,  terms  of  E are 
r-5  — 


4ire  r 


- (r  -2cos  0 + 0 sin  0)  U(t  - tQ ),^  (A. 19) 

3 — — c 


and  this  is  exactly  the  static  field  of  a dipole  q;  the  field  arises  at 
t = t and  propagates  outward  with  velocity  c.  This  is  the  change  in 
the  static  field  due  to  the  current  JQ,  as  noted  above. 

The  above  fields  have  been  derived  for  an  electric  current  im- 
pulse, but  if  electric  and  magnetic  quantities  are  properly  interchanged 

they  become  the  fields  of  a magnetic  impulse.  Discussions  of  magnetic 

* 

sources  are  found  in  the  literature,  and  the  application  to  the  impulse 
current  follows  exactly  the  lines  indicated  above. 

As  an  example  of  the  use  of  these  results,  let  us  find  the  field  of 
an  oscillating  dipole  which  is  suddenly  shut  off  at  t = t^.  We  shall 
confine  our  remarks  to  Eq  because  that  component  has  all  the  essential 
characteristics.  The  current  is  represented  by 


J(t)  = z q cos  co  t S (Q)  U(t^  - t). 


(A. 20) 


The  electric  field  is  given  by  eq.  (2.6)  of  the  text: 


* See,  forexample,  Schelkunoff,  S.  A.,  Electromagnetic  Waves,  chap.  VI. 
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00 

:e(Qt)  = / / J(PtQ)  G0(Qt,  PtQ,  z)  dtQ  dvp.  (A. 2 1 ) 


V -03 


The  0-component  of  GQ  is  given  by  eq.  (A. 17).  This  gives 


E0  (Qt)  = 


q sin  6 n I 


:os^t0  U(t,  -tQ)  ffD(t-t0-|) 


+ 72  8 (*  " *o  -7)  + 73  U(t  - t„  - f)  dtD. 


(A. 22) 


The  first  integral  on  the  right-hand  side  is  integrated  by  parts: 


[cos  63 10  U(t1  - tQ)]  D(t  - t0  - J)  dto 


- COS  63  tQ  U(tj  - tQ)  S (t  - tQ  - J) 


tQ=  -00 


3 sin  a;  tQ  U(tj  - tQ)  + cos  63  tQ  S (tj  - tQ)  5 (t  - tQ  --)  dtQ 


= 0 - 


sin  63  (t  - -|)  U [(tj  + f)  “ t]  - coswtj  S [t  - (tj  + -)] 


(A.23 ) 

The  first  term  is  zero  because  the  8 -function  is  zero  at  ± 00  . The 
second  term  is  the  "radiation”  field  of  the  dipole.  It  is  sinusoidal  in 
time  and  ends  abruptly  at  (t^  + E).  The  third  term  is  a pulse  which 
radiates  at  the  instant  the  current  is  stopped.  It  has  a magnitude  pro- 
portional to  cos  63  tj ; i.  e.,  to  the  instantaneous  value  of  the  current  when 
it  was  stopped.  If  the  current  were  zero  at  tj  there  would  be  no  pulse, 
because  no  infinite  acceleration  would  be  needed  to  bring  the  charge  to 
rest. 
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The  second  integral  in  Eg,  eq.  (A. 22),  is  integrated  directly. 


J COS  Cl 1 tQ  U(tj  - tQ)  3 (t  - tQ  - — ) dtQ 
-oo 

= cos  o(t  - f)  U [(tj  +f)  - t]  . 


(A.  24) 


This  term  (when  multiplied  by  its  factors)  is  the  -L  part,  the  induction 

rc. 

field,  of  the  dipole.  It  also  ceases  abruptly,  but  there  is  no  pulse  when 


it  stops. 


The  third  integral  in  Eg  is  integrated  by  parts. 


La>cosa)to  U(*l  - to)  U(*  - *o  -f)  dto 
= 7J  sin  " tQ  U(tj  - tQ)  U(t  - tQ  - -r) 


*0  = -°° 


+ f00^sinwto[U(tl  -to)S(t-t0-f)+  8(tj  -tQ)  u(t  - tQ  - I)] 


sin co  tr 


U(tj  -t0)U(t-t0-|)  +^sin<w(t-|)u[(t1  + £)-t] 


t = —00 

xo 


+ ~ sin  co  tj  u[t  - (tx  + -|)] 


(A. 25) 


This  second  term  is  the  customary  i portion  of  the  field  o^L  an  oscil- 

r j 

lating  dipole,  and  it  stops  at  t = (ti  +— ).  The  last  term  is  the  residual 

c 

static  dipole  field  existing  for  t > (tj  + -^-).  It  is  just  a constant  extension 
of  the  time  dependent  part;  i.  e.,  the  second  term,  which  is  cut  off  at 
t = (t,  +— ).  The  first  term  is  zero  at  the  upper  limit  and  at  the  lower 
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limit  may  be  interpreted  as  the  static  dipole  field  which  existed  before 


the  current  was  started.  (The  transients  which  occurred  when  the  cur- 
rent was  started  have  been  ignored  because  it  may  be  assumed  that  they 
have  all  died  down  by  time  tj.)  Both  of  these  static  fields  may  be  ignored, 
along  with  the  static  portion  of  $ mentioned  above. 

Thus  we  have,  for  the  nonstatic  portion  of  the  0-component  of 
electric  field  due  to  an  oscillating  dipole  which  cuts  off  at  t = tj : 


co 


= 3-SiSJ  IE 

4it  V e L rc 


sinw 


(t  - £)  + _L  cos  cu  (t  - -Jr)  + — \ sin  co( t - ~) 


c'  2 
r 


COT 


{U  [(ti  + f ) - t|  - q ~ cos  « tj  8 [t  - (tj  + £)] 


(A. 26 ) 
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APPENDIX  II.  CAVITY  GREEN«S  FUNCTIONS 


Z: 


1.  INTRODUCTION 

This  appendix  is  concerned  with  the  expression  of  electric  and 
magnetic  fields  in  a closed  region  in  terms  of  their  values  on  the 
surface.  The  fields  will  be  written  as  a linear  combination  of  nor- 
mal modes  appropriate  to  the  region;  these  modes  satisfy  the  same 
differential  equation  as  the  field,  and  certain  boundary  conditions. 

The  result  is  a generalized  Fourier  series,  where  the  function  to  be 
expanded  is  known  only  on  the  surface. 

This  problem  has  been  discussed  in  the  literature,  for  both 
scalar  and  vector  fields.23’  3010  33  It  has  been  shown34  that  the  modes 
obtained  from  a cavity  with  a homogenous  "short-circuit"  boundary 
condition  (i.e.  , n xE  = 0 onS)  form  a complete  set  for  the  expansion 
of  vector  fields  with  vanishing  tangential  component  on  the  surface. 
The  earlier  writers  on  the  present  subject  assumed  that  these  short- 
circuit  modes  were  complete  for  the  expansion  of  fields  with  nonzero 
tangential  component  on  S,  but  it  turns  out  that  they  are  not,  and  ad- 
ditional terms  must  be  introduced.  This  lack  of  completeness  was 
first  pointed  out  by  Wigner.*  It  was  independently  discovered  by  the 
author  when  an  attempted  expansion  failed  to  yield  the  correct  results. 

* Private  communication.  Professor  Wigner  and  Dr.  Teichmann  very  kindly  made  some  as  yet 
unpublished  material  available  to  the  author. 
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The  function  Teichmann  and  Wigner  use  to  complete  the  expansion 
is  different  from  those  used  here,  although  it  will  be  shown  that  the 
two  methods  are  equivalent.  The  complete  expansions  presented 
here  were  originally  developed  by  Schwinger,  although  the  present 
method  of  derivation  is  somewhat  different.  In  the  same  paper  with 
the  complete  expansions  Schwinger  produced  an  alternative  (and  ap- 

3 1 

parently  simpler)  expansion  using  only  the  short  circuit  modes. 

2.  EXPANSIONS  IN  EIGENFUNCTIONS 


Let  V be  a homogeneous  region  bounded  by  S,  with  outward 


normal  n.  There  is  an  external  source  with  a time  dependence 
producing  a field  in  V,  and  the  electric  and  magnetic  fields  in  V sat- 
isfy the  wave  equation 

VxVxE-k2E  = 0.  (A.  27) 

A knowledge  of  n x E or  rt  x H on  S determines  E and  H in  V. 

Now  define  the  functions  Ea  as  solutions  of  the  equation 

Vx  VxEa  - l^2  Ea  = 0 (A.  28) 

subject  to  the  boundary  condition 

nx  Ea  = 0 on  S.  (A.  29) 

To  show  the  orthogonality  of  these  functions,  use  the  vector  analogue 
of  Green’s  theorem:36 

j (A'?xVx  B - B*?x7x  A)dv 

J 


V 


= J (BxVxA-AxVx  B)  ■ n ds , 


(A.  30) 
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X 


If  A = E-  and  B = Eg  , then  from  eqs.  (A.  28)  and  (A.  2 9), 

(ka2-kpZ)f  Ea.Egdv=0.  (A.  31) 

J V “ ~ 

The  functions  corresponding  to  different  eigenvalues  are  therefore 
orthogonal,  and,  if  there  is  degeneracy,  appropriate  linear  combina- 
tions will  yield  an  orthogonal  set  of  functions  going  with  the  same 
eigenvalue.  The  functions  may  also  be  normalized,  by  setting 

f Etx  ' Eg  dv  = . (A.  32) 

J y~  ~ 

We  also  need  the  magnetic-type  modes,  denoted  by  1^  and 
defined  by 

ka  Ha  = V x Ea  . (A.  33) 

Then,  from  eq.  (A.  28), 

ka  Ea  = . (A.  34) 


The  functions  Ha  satisfy  the  wave  equation 


VxVx^  - C Ha  * 0 


(A.  35) 


with 


n x V x H =0  on  S. 
a 


(A.  36) 

To  demonstrate  their  orthogonality  use  the  identity 

V • (A  x B)  = B • V x A - A • \7  x B . (A.  37) 

Thus , 

/ H*  * Hg  dv  = —f  Ha  • V x Eg  dv 


H 

V —P  kg*  V 


=—  \(*'  (^x!a)dv+  /js" dv'l 

k^V  V (A.38)J 
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Substitute  (A.  34)  and  use  the  boundary  conditions  to  obtain 


f Ha  • H 

J 


'/3 


dv  = 8 


a/3 


(A.  39) 


The  sets  of  functions  Ea  and  Ha  are  the  electric  and  magnetic 
fields  - the  normal  modes  - which  can  exist  in  the  closed  cavity. 
Note  that,  because  of  the  normalization,  Ea  and  Ha  do  not  form  an 
electromagnetic  field;  if  Ea  is  an  electric  field  then  J "X  Ha^  is 
the  associated  magnetic  field. 

Now  consider  the  eigenvalue  zero.  It  must  be  examined  sepa- 
rately because  the  above  definitions  are  not  valid  in  this  case.  Let 
the  index  p denote  functions  going  with  this  eigenvalue,  and  from  now 
on  reserve  a for  eigenvalues  other  than  zero.  The  electric-type 
functions  satisfy 


with 


V x VxEp  = 0 


n x E„  = 0 on  S. 
_ _ P 


(A.  40) 


(A. 41) 


From  these  two  equations  it  is  possible  to  show  that  E^  is  irrotational. 


Thus,  set 


Then 


V x Eq  = Vu 


(A.  42) 


V u = 0 


(A.  43) 
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The  only  solution  to  (A. 43)  and  (A.  44)  is  u = constant,  and  hence 

V x Ep  = 0.  (A.  45) 

The  functions  En  are  irrotational,  and  all  the  functions  E are 

r 

solenoidal.  Inasmuch  as  the  fields  to  be  expanded  are  solenoidal,  it 
might  seem  reasonable  to  assume  that  the  functions  E^  are  unneces- 
sary here,  although  they  may  be  useful  for  the  expansion  of  a static 
field.  However,  to  obtain  a complete  set  it  is  necessary  to  have  the 
eigenfunctions  going  with  all  eigenvalues;  therefore  we  shall  keep 
the  irrotational  functions. 


There  is  an  infinite  number  of  functions  5>  • They  may  be 
defined  by 

Ep  = — V^p  , (A. 46) 

v 

where 

V%p  + \p%p  = 0 in  V (A. 47) 

and 

<j>p  , « 0 on  S.  (A. 48) 

By  using  Green's  theorem  it  is  shown  that  the  functions  </>p  form  an 
orthonormal  set: 

J <t>p  <PV  dv  = Spv  » (A. 49) 


and  it  then  follows  that 


f Ep  ‘ Ej, 

* 


dv  = S 


pv 


The  condition  n x Ep  = 0 on  S is  also  satisfied. 

Similarly,  the  functions  H Q are  defined  by 

Up  Wp  , 

XP 

where 

+ \p\  = 0 in  V 


(A.  50) 


(A. 51) 


(A.  52) 


The 


^ p 

= 0 on  S . 

3 n 


Kp's  may  be  normalized: 

H 


and  then 


(A.  53) 


(A.  54) 


(A.  55) 


Also,  n - H = 0 on  S.  This  is  the  appropriate  magnetic-type  bound- 
ary condition  to  satisfy  here,  rather  than  (A. 36),  because  V * is 
identically  zero. 

The  functions  <pp  and  i pp  are  complete  sets  for  the  expansion  of 
scalar  functions  (which  satisfy  certain  requirements  as  to  continuity, 
etc).  They  are  the  functions  used  by  Sommerfeld  in  his  discussion 

30 

of  the  scalar  cavity  Green's  function.  From  the  theory  of  Laplace's 
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equation  it  can  be  shown  that  there  are  no  functions  going  with  = 0, 
(except  0 = constant).  Therefore,  none  of  the  functions  E^  and  H ^ 
are  possible  source-free  static  fields,  because  their  divergence  is 
not  zero. 

Now  assume  that,  at  every  point  within  the  cavity,  it  is  possible 
to  write  E and  H as  linear  combinations  of  the  above  functions: 


E = 2 aa  Ea  + 2 a.p  E p 
— a — P H K 

(A.  56) 

H = J Ha  + 2 bp  H,  . 

(A.  57) 

By  the  orthogonality  relations, 

a = f E • E 
a I _ _a 

* * \r 


dv 


xp-f  E ■ 5°  dv> 


(A.  58) 


and  similarly  for  the  b's.  The  expansion  coefficients  are  found  in 

terms  of  n x E from  Green's  theorem.  From  (A.  27)  and  (A.  28)., 
2 

- k ) 

V "V 


(ka2  - k2)  f E • Ea  dv  = f (E  • V x V x £a  - Ea  • V x V x E)  dv 


= f (Ea  x V x E - E x V x Ea)  • n ds. 
JS  ~ — ~ - 

(A..  59) 


Therefore , 


/„£  • 5“ dv  * L (n  x E)  • Ha  ds.  (A.  60) 


k2-ka2  -s 


t 
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Similarly 


/ ^ 
4 V 


E 0 dv  = 0 


f H H dv  = — / (n  x E)  • Ha  ds 

J V~  ~ k2  - k 2 •'S 


/ H.  Hpdv  = — f (nxE)  .H, 
*V  OJ  fj.  J S 


ds. 


Therefore  , 


(A. 61) 
(A.  62) 

(A.  63) 


E = Y — Ea  f (n  x E)  • Ha  ds  (A. 64) 

- V k2  - ka2  - •'s 

H =V  . Ha  f (nxE).H  ds+V-^f(nXE)  ■ H^ds. 

“ Vk  - ka  JS~  - -rjs-  - -r 

(A.  65) 


The  solenoidal  set  Ea  alone  is  sufficient  to  represent  E but 
both  sets,  Ha  and  H 0,  are  needed  for  H.  This  difference  is  due  to 
the  different  boundary  conditions  on  the  Ea  and  Ha.  The  expansion 
is  analogous  to  a Fourier  series  over  the  half-interval  0<  x<it.  The 
sines  alone  form  a complete  set,  but  the  cosines  must  also  contain 
the  constant  term.  The  constant  term  is  derived  from  the  eigenvalue 
zero. 
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32;, 33  . 

Teichmann  and  Wigner  complete  the  expansion  by  adding 

an  ir rotational  vector,  as  follows* 


E = 2 a E (A.  66) 

— a — 

H = I ba  Ha  + Vu,  (A.  67) 

— a — 

where 

= 0 in  V (A.  68) 


= n • Hr  on  S.  (A.  69) 

d n 

The  function  u is  defined  by  eqs.  (A.  68)  and  (A. 69),  in  terms  of 
normal  H on  the  surface.  Now  the  coefficients  aa  and  ba  are  given 
in  terms  of  tangential  E,  and  it  would  seem  advantageous  to  com- 
pletely specify  the  expansion  in  terms  of  tangential  E,  as  was  done 
above,  in  eqs.  (A.  64)  and  (A.  65).  Actually,  of  course,  tangential  E 
specifies  normal  H,  so  there  is  no  essential  simplification.  The  two 
methods  are  in  fact  equivalent,  as  will  now  be  shown. 


written  as 


30 


N(Q,P)  t p (Q)  *p{ P)  • 

p V 


(A. 71) 


When  this  expansion  for  N(Q,P)  is  substituted  into  (A.  70)  the  gradi- 
ent may  be  taken,  and  it  may  then  be  shown  that  the  result  is  equal  to 
the  second  series  in  eq.  (A.  65). 

The  procedure  outlined  above  does  not  work  when  k = ka_;  i.  e.  , 
when  the  actual  field  oscillates  at  a resonance  frequency  of  the  cavity. 

In  such  a case,  since  f E • E dv  is  finite,  the  surface  integral 

JV~  ~ 

f n x E • Ha  ds  must  be  zero.  This  is  the  exceptional  case  when  the 

J s-  — - 

uniqueness  theorem  is  not  valid.  A more  precise  uniqueness  theorem 
would  read:  values  of  nx  E on  a closed  surface  uniquely  determine 
the  (source-free)  interior  field  unless  the  oscillation  frequency  hap- 
pens to  coincide  with  a resonance  frequency.  In  this  case,  since  the 
resonant  mode  has  n * E = 0 on  S , it  may  be  added  to  the  field  with- 
out changing  n x E on  S. 

Since  the  functions  Ea  correctly  represent  E they  must  also 
be  adequate  to  represent  H.  We  may  therefore  write 


H =J  ca  Ea 
— a . — 


where 


C(X  2 


f (n  x H)  . H ds. 

_ u c J c — 


(A.  72) 


(A.73) 


k - ka“  " S 

H is  now  expressed  as  a sum  of  solenoidal  functions,  in  terms  of 
n x H on  S.  However,  if  n x H (rather  thanuxE)  is  known  on  S, 
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the  expansion  for  E would  use  both  sets  of  functions,  and  Hp.  In 
the  present  case  the  theory  is  applied  to  a radiating  cavity  where 
n x E is  known  (or  assumed)  in  the  aperture;  therefore  the  first  ex- 
pansions above  are  used.  Other  combinations  would  be  possible  if  we 
had  followed  Slater  and  allowed  a piece  of  the  surface  to  be  an  open 
circuit  (magnetic  conductor).  If  it  covered  the  entire  surface  then 
the  Ha 's  would  be  the  Ea's  of  the  above  discussion,  and  obviously 
nothing  new  would  be  gained.  If  mixed  boundary  conditions  were  al- 
lowed then  both  sets  of  functions  would  be  needed,  in  general. 

It  is  of  interest  to  show  that  the  series  (A.  64)  and  (A.  65)  satisfy 
the  Maxwell  equations.  However,  as  pointed  out  by  Slater,  the  series 
for  E and  H cannot  be  differentiated  term-wise;  but  it  can  still  be 
shown,  for  example,  that  the  expansion  for  curl  E agrees  with  the 
expansion  for  H.  Curl  E may  be  expanded  as  a sum  of  the  Ha  and 


H p functions: 


V x E = S Ha  f Ha  • V x Edv  + IH  f Hp  • V x E dv. 
a J , , p J it 


(A.  74) 


Use  the  divergence  theorem  to  obtain 


V * 2 S E • Ea  dv  + j (n  X E)  • Ha  ds 


l%,f  <?*  B -2 p ds- 


(A.  75) 
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The  first  part  of  the  first  sum  in  (A.  75)  is  equivalent  to  what  one 
would  obtain  by  taking  the  termwise  curl  of  (A.  64).  Substitution  of 


(A.60)  yields 


V x E = - jc Ofj.  2 H,  — 
— a — , c 


k - k„ 


f (nxE)  . Ha  ds 

J C* 


j r 1 

2 Hp J (nx  E)  • Hp  ds 

^ co/xS 


(A.  76) 


The  bracketed  expression  is  exactly  the  expansion  for  H,  eq.  (A.  65). 


2 . 

The  cavity  Green's  function,  r(Q,P,k)  is  defined  by 


H(Q)  = f r(Q,P,k)  K (P)  ds 

« A 


(A.  77) 


A comparison  of  this  with  eq.  (A.  65)  shows  that  ^ may  be  written  as 

r(Q,P,k)  =y  ■~U€— - Ha  (Q)  Ha  (P)  +Y  j-H  (Q)H  (P).  (A.  78) 
2 , 2 — — L scou  H H 


bbV  - V ^ ~ 

CL  K-  p 

This  is  the  tensor  (dyadic)  cavity  Green's  function.  Note  that 

2 T 

r2(Q.p)  = r (p , q) . 

3.  APPLICATION  TO  WAVEGUIDE  CAVITY 


(A.  79) 


These  results  will  now  be  applied  to  a cavity  in  the  shape  of  a 
waveguide.  The  cross  section  is  arbitrary  but  uniform  along  the 
z-axis;  there  is  a short  circuit  at  z = 0 and  the  aperture  is  at  z = ^ 
(Fig.  3 9). 
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a.  CONSTRUCTION  OF 
EIGENFUNCTIONS 


Fig.  39.  Waveguide  cavity. 


The  first  thing  to  do  is  to  set 
up  the  cavity  eigenfunctions.  The 
solenoidal  functions  may  be  con- 
structed by  observing  that  they 
consist  of  two  sets,  one  TE,  and 
one  TM,  to  the  z-axis.  The  cross- 
sectional  variation  will  be  just  the 
waveguide  modes,  described  in 
section  II,  and  the  z-dependence  is 


p 

such  as  to  make  E zero  at  z = 0 and  z = The  TE  modes  may 


therefore  be  written  as 


kJE  (xyz)  = Nav  e„TE  (xy)  sin 


a = 1,  2,  3, 


(A.  80) 


Up  to  now  there  has  been  a single  subscript,  a , on  the  mode  functions. 
This  actually  stood  for  three  indices,  because  the  cavity  is  a three- 
dimensional  region.  In  eq.  (A.  80),  a takes  on  integral  values,  and  the 
subscript  v stands  for  two  indices. 

The  normalization  condition  determines  Na,v: 


2 S'  j aiTZ  r 

- f'W  f sin  — — dz  J ej7E  (xy)  • evTE  (xy)  ds.  (A.  81) 

* n * A “ “* 


t 
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dirz 


With  the  result  that 


% vE  fcyz)  ~J~  sin 


l ~ 


TE  / v 
ev  (xy)  , 


a = hi  2, 3.V-  (A.  82) 


TE 

(xyz)  = 0. 

Then,  from  eq.  (A.  33), 


a = 0. 


(A.  83) 


V x £ 


(A.  84) 


air  airz 


/ cos (-xev  +ye  )+  z component.  (A.  85) 


The  main  interest  here  will  be  on  the  transverse  fields,  so  the  z- 
component  is  merely  indicated.  The  quantity  in  brackets  is  equal  to 
z x e , and,  from  eq.  (2.46),  the  TE  magnetic  modes  become 


H av  (xyz)  = / cos hj ^ (xy)  + z cpt ,.  a=l,2,3,***  (A.  86) 

“ 'Uav  l ~ 


HoJ E (xyz)  = °» 


a = 0. 


(A.  87) 


Similarly, 


nr-  a it  z 

EttJM(xyz)  =J—  cos — — hJM(xy), 


a = 1,2,3,  (A.  88) 


^ h™  (xy), 


a = 0 


(A.  89) 
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Now  (A.  93)  and  (A.  94)  are  just  the  equations  satisfied  by  the 

8 6 

potential  functions  which  generate  the  TE  waveguide  modes.  Hence 
the  eigenvalues  kv  must  be  the  cutoff  wavenumbers  for  the  TE  modes, 
and  the  functions  Ty  (xy)  are  equal  to  the  potential  functions,  apart  from 
a constant.  This  furnishes  an  opportunity  to  express  the  functions 
Hq  in  terms  of  . 

The  potential  functions  F^  are  defined  by 

sTE(xy)  = ^2 x 5/  fry)*  9?) 

and  they  have  a z-component  only.  If  the  curl  is  expanded  it  is  seen 
that 

STE(xy)  = -z  x V2Fv  (xy),  (A.  98) 

where  Fy  is  the  z-component  of  F^.  If  we  substitute  for  e^E  from 
eq.  (2.46),  we  have 

hvTE(xy)  = V2Fv(xy).  (A.  99) 

To  find  Ty  in  terms  of  F , set 

Tv  (xy)  = Nv  Fv  (xy).  (A.  100) 


\ 

V - 
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Then 


lmf*£  dv  = /’ 


J , P*z 


cos  

0 t 


dz 


/ Tv  (xy) 

» A 


ds 


l N* 


-t  -if. 


2 k;  'A 


FvVZFv  ds 


* 2 1 


= ,_N2—  f hTE. 

O V ..2  J A -v 


«l  ” K 


hTE  ds. 


Therefore , 


NV  = ,/— 


and 


= 0. 


Then 


Hpv  = 


\p 


v 


i r ^itz  i 

V'/'pv= v NVP^  (xy)  cos— — 

Xpv  L ^ J 


PTTZ 

T 


4- 


pTTZ 


COS h^E (xy)  + z^cpt. , p-  1 , 2 , 3 , • • • (A.  105) 


(A.  101) 


p - 1, 2,3,  •”  (A.  102) 


(A.  103) 


(A.  104) 


Xpy  ^ 


Hoj/(xyz) 


■£ 


i - 


h JE  (xy) 


P=  0. 


(A.  106) 


The  functions  in  eqs.  (A.  105)  and  (A.  106)  are  the  irrotational  functions 
needed. 
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b.  EXPANSION  OF  MAGNETIC  FIELD 


In  section  III,  in  connection  with  the  normal  modes  of  the  radiat- 
ing waveguide-cavity,  it  was  stated  that  the  transverse  magnetic  field 
obtained  from  the  eigenfunction  expansion  agrees  with  that  computed 
from  transmission  line  theory.  This  will  now  be  shown. 

The  expansion  is  given  in  eq.  (A.  65),  and  now  replace  cj  by  the 
complex  frequency  p.  The  integrands  are  zero  except  in  the  aperture, 


so 


jpe 


H=HH„ 


" **-aJV-k2  "A 

av 


[ 
J A 


PH  A 


f (zx  E) 

J A 


H ds. 
-Ipv 


(A.  107) 


Consider  first  the  TM  modes  of  the  first  sum.  Substitute  for  H 


TM 


av 


from  eqs.  (A.  88)  and  (A.  89),  retaining  only  the  transverse  fields. 

(-l)tt  f (z  x E)  »hJMds 

• A 


t m TM  2 a.TTz  jpe 

H (xyz)  =2  2 h'y^(xy)  — cos 

a=l  v l l kZ  _ k2 


av 


1 jpe 

.r,TM/  X 

+ 2 \ (xy) T 

v ik2_k2 


f (z  * E)  • h™  ds.  (A.  108) 

* A 


ov 


a 


The  (-1)  comes  from  the  cos  SJIZl  in  the  integrand,  evaluated  at 

& 


z = 't  . The  notation  is  rather  cumbersome,  so  the  following  list  of 


symbols  is  included  for  clarification: 


a = number  of in  the  waveguide  (for  the  eigenfunctions  Ea,  etc.) 

2 


v = waveguide  mode  number  (stands  for  two  indices) 
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T? 


kv~  Zv  x cut-off  wavenumber  of  v±h  waveguide  mode 
k = pjjle  =2rir  fjle  x complex  oscillation  frequency  of  the 


normal  mode  of  the  radiation  problem, 
k =--  2ir  JjjT  x oscillation  frequency  of  the  avth  cavity  mode. 


+ k,.  mode. 


r 2 = - k2  • 


Hence, 


k2  - kiv  = - y}  -(-J-) 


(A,  109) 


(A.  110) 


Rearrange  eq.  (A.  108)  to  give 


H(xyz)  — h™  (xy)  f E ■ e™  d,(-L 


r- ' " a 

At1'11 

i-i  U7 


X cos  

2 l 


(A.  Ill) 


But,  as  may  be  verified,  the  quantity  in  the  brackets  is  the  cosine  ex- 

, / coshT™  z \ „ 

pansion  of  ( ) , over  the  range  0 ^ z $ v.  Also,  for  TM 

' sinh  yv  l ' 

modes,  the  quantity  is  the  waveguide  characteristic  admittance, 

y-^  ^ (see  section  II).  Hence  the  contribution  of  the  TM  modes  to  the 
transverse  magnetic  field  reduces  to 

, cosh7j,  z „ 

H (xyz)  = 2 -yv™  h v™  (xy)  / E • 5,™  ds.  (A.  112) 

v sinhyv  -i  A 


Now  consider  the  TE  modes.  We  will  consider  the  irrotational 
modes  together  with  these  because  their  transverse  dependence  is  of 
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*he  TE  type.  Substitute  for  H^p  and  Hpv  in  (A.  107)  -to-obtain 

~ 2 


H (xyz)  =I£  h TE  (xy) 
a = l v ~v 


" air  " 

2 

l 

l 

\ 

airz 
cos  (- 1) 


a JP^ 


k2  - k2  " A 

av 


f (f*E)-hJEds 


00 


tSS  V™liir)  - 

Pm  1 ^ ^ 


2r*v  1 


piTz 


cos (-if  — [ (z*E)’  hJ 

4 » » 


TE 


ds 


PM  A 


+ 2 hTE(xy)  -—  [ (z*E)  . hTE  ds. 
v ~v  / PM  JA  - — -v 


(A.  113) 


After  using  eqs.  (A.  96),  (A.  109)  and  (A.  110),  the  equation  becomes 


H (xyz)  = 2 hj* E (xy)  f E •'  e 

V J — 


TE  I 1 J » ±j_  (-1)P 


pirz 


CQS- 


{Pf  "-up"  «lijY  * 


30 


+cr4viPe(~1> 


a (¥)‘ 


(-1) 


afrz 

cos y.  (A.  114) 


The  last  two  sums  may  be  combined  (change  the  dummy  index  p to  a ) 
and,  after  a little  manipulation,  one  obtains 


t ^ u 

H (xyz)  =2  hvTE(xy)  f E • 

JPM  A 


TE  , 
% ds- 


00 


7 


V 


+ 2 - (-if 

a=l  £ , ,air 


O-TTZ 


COS 


7 


(A.  115) 
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The  bracketed  quantity  is  again  ( —r)  5 and  is  the  charac- 

''sinh7v/v'  'jfiuV 

teristic  admittance  for  the  vth  TE  mode.  Therefore,  the  contribution 
of  the  TE  and  irrotational  eigenfunctions  has  the  same  form  as  that 
of  the  TM  eigenfunctions.  Let  the  subscript  v denote  all  the  wave- 
guide modes,  and  as  a final  result,  write 


t cosh yv  z . 

H(xyz)  =Z~yv  h^xy)  / E ’ ez 

v sinh  Jv J A 


(A.  116) 


Now  consider  the  cavity  from  waveguide  theory.  We  have  a sec- 
tion of  waveguide,  shorted  at  z = 0,  and  with  z x E prescribed  at 


z = % . F rom  section  II, 

H (xyz)  = £ I (z)  h (xy) 

— v " — " 

E(xyz)  * 2 V (z)  e (xy)  . 

— v — 

At  any  point  z,  lv  is  the  sum  of  two  waves: 

I„  W - Av  (e**  + 

Let  z = -C,  then 


(A.  117) 
(A.  118) 


(A.  119) 


(A.  120) 


M*> 

A v-  (A.  120) 

2 cosh  yv 

But,  from  eq.  (3.27), 

lvW  = Yv  coth  yvl  - - yv  coth  yvl  J E • evds.  (A.  121) 
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It  follows  then  that 


coshy^z 

Mz)  = -vv  [ E‘iv 

— T™  0 " A 


sinhyv^  A 


t cosh  yv  z r 

H(xyz)  = 2 -yv  (xy)  I E • e v ds. 

v sinh-yv^  A 


It  has  therefore  been  shown  that  the  transverse  magnetic  field 


of  a waveguide  cavity  calculated  from  the  expansion  in  cavity  eigen- 


functions agrees  with  the  much  simpler  calculation  based  on  waveguide 


theory. 
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APPENDIX  III.  THE  STATIONARY  PRINCIPLE 


It  is  desired  to  show  that  the  propagation  constant  k,  as  calcu- 
lated from  eq.  (3.17)  of  the  text,  is  stationary  with  respect  to  first- 
order  variations  of  tangential  E about  the  true  value.  This  stationary 
principle  follows  from  Schwinger's  work.  This  particular  proof  is 
due  to  Professor  Rumsey  (Final  Engineering  Report  400-11,  ARDC). 

The  equation  under  investigation  is 


(A. 124) 
(A. 125) 


Take  the  first  variation: 


f j S K(Q)  T(Q,P,k)K  (P)  dsp  dsQ 


AQ  ap 


+ f f K(Q)  Sr  (Q,P,k)K(P)dspdsQ 

•'Aq  ;Ap 

+ J J K(Q)r  (Q,p,k)s  K:(F)  dsp  dsQ  = 0.  • (A. 126) 


aq'ap 


1 2 

Now  r and  T are  symmetric,  and  therefore  T is  symmetric  (with 


interchange  of  Q and  P).  The  last  term  therefore  becomes  (c.  f. 
eq.  (3.14)  of  the  text) 
j j K(Q)  T (Q,P,k)  S K (P)  dsp  dsQ 


Aq  Ap 


= J J S K(P)r(P,Q,k)K  (Q)  dspdsQ.  (A.127) 


Aq  ¥ Ap 
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If  the  order  of  integration  is  interchanged  this  is  just  equal  to  the  first 
term  of  (A. 126).  But  the  first  term  is  zero  when  evaluated  at  the 
correct  K,  because 

z x f r (Q,P,k)K  (P)  dsp  = z X A(Q)-H(Q)  =0.  (A.134) 

" A L 


Tangential  H is  continuous  across  the  aperture. 


Thus  the  first  variation  is 


f f K(Q)S  T (Q,P,k)K  (P)  dsp  dsQ  = 0,  (A. 135) 

•I  A _ J A 


Aq  j Ap 


8k  f f S(Q)5L  K(P)  dsp  dsQ  = 0.  (A. 136) 

J A J A 3 k 


AQ  AP 


The  integral  is  presumably  not  zero,  so 


8 k = 0. 


(A. 137) 
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